AD-A181  17S 


TGAL-86-5 


MODELS  OF  THE  FREQUENCY  DEPENDENCE 
OF  Q  IN  THE  MANTLE  UNDERLYING 
TECTONIC  AREAS  OF  NORTH  AMERICA, 
EURASIA  AND  EASTERN  PACIFIC 


Z.A.  Der,  W.W.  Chan,  A.C.  Lees,  and  M.E.  Marshall 
Teledyne  Geotech  Alexandria  Laboratories 
314  Montgomery  Street 
Alexandria,  Virginia  22314-1581 


OCTOBER  1986 


prrds 

electe 

JUN  o  8  198/1 


V 


FINAL  TECHNICAL  REPORT: 
ARPA  ORDER  NO:  4511 


PROJECT  TITLE:  Frequency  Dependence  of  Q 
CONTRACT:  F08606-85-C-0023 


Approved  for  Public  Release;  Distribution  Unlimited. 


Prepared  for: 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


Monitored  By: 

AFTAC/TGR 
PATRICK  AFB 
FLORIDA  32925-6001 

The  views  and  conclusions  contained  in  this  report  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  representing  the  official  policies,  either  expressed  or  implied,  of  the  Defense  Advanced 
Research  Projects  Agency  or  the  U.S.  Government. 


,w mw 


UNCLASSIFIED  _ 

SECURITY  CLASSIFICATION  OP  THIS  PAG~ 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


Form  Approved 
OMBNo.  0704  01 88 
Exp  Date  Jun  30.  1986 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S)  ' 

TGAL-86-5 

6a  NAME  OF  PERFORMING  ORGANIZATION 

6b  OFFICE  SYMBOL 

Teledyne  Geotech 

(If  applicable) 

Alexandria  Laboratories 

6c.  ADDRESS  (City.  State,  and  ZIPCode) 

314  Montgomery  Street 

Alexandria,  Virginia  22314-1581 

8a.  NAME  OF  FUNDING /SPONSORING 

8b  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

DARPA 

GSD 

8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

1400  Wilson  Blvd. 
Arlington,  VA  22209 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited . 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFTAC/TGR 


7b  ADDRESS  (City,  Slate,  and  ZIPCode) 

Patrick  Air  Force  Base 
Florida  32925-6001 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F08606-85-C-0023 


10.  SOURCE  OF  FUNDII)lG  NUMBERS 


PROGRAM 
ELEMENT  NO 

62714E 


PROJECT 

NO 


I  WORK  UNIT 
ACCESSION  NO 


PA  DT/ 5122 


11.  TITLE  (Include  Security  Classification) 

Models  of  the  Frequency  Dependence  of  Q  in  the  Mantle  Underlying  Tectonic  Areas  of  North 
America,  Eurasia  and  Eastern  Pacific 


12.  PERSONAL  AUTHOR(S) 

Z.A.  Der,  W.W.  Chan,  A.C.  Lees  and  M.E.  Marshall  _ _ 


13a  TYPE  OF  REPORT  1 13b  TIME  COVERED  ^4  DATE  OF  REPORT  (Year.  Month.  Day)  15  PAGE  COUNT 

Final  Technical  Report  from  Mar  85  TO  Sep  86  1986  October  30 _  171 


16.  SUPPLEMENTARY  NOTATION 


17  COSATI  CODES  18  SUBJECT  TERMS  (Continue 

FIELD  GROUP  SUB-GROUP 

- _ - - - Attenuation,  Freq> 

- 2® - - Eurasia,  North  Am' 

19.  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


Attenuation,  Frequency  Dependence  Q,  Magnitude, 
Eurasia,  North  America,  E.  Pacific 


~  The  results  of  wide-band  data  analyses  for  the  determination  of  frequency  dependent  Q  models 
of  the  mantle  under  the  tectonic  areas  of  Eurasia,  North  America  and  the  Eastern  Pacific  are  presented. 
For  tectonic  areas  of  Eurasia,  the  spectra  of  short  period  P  waves  from  Kazakh  explosions_are  used  to 
estimate  t*  in  the  1  to  8  Hz  range.  (For  shield  to  tectonic  paths,  tp  -  0.35  sec  which  gives  tp  -  0.6  sec 
for  tectonic  to  tectonic  paths.  P  witve  rise  times  from  the  same  events  give  upper  bounds  on  the  esti¬ 
mates  of  to  be  -  0.8  sec.  The  large  difference  between  the  spectral  and  rise  time  estimates  may  be 
due  to  scattering.  In  the  0.5  to  1.3  Hz  range,  PP  and  P  amplitude  ratios  give  -  0.5  sec  for  PP’s 
relative  to  surface  reflection  points  under  the  Eurasian  shield^  Amplitude  ratios  of  long  period  multiple 
ScS  phases  were  used  to  estimate  attenuation  in  the  0.3  to  0.9  Hz  band.  Under  tectonic  Eurasia,  these 
analyses  gave  ts*cS  ~  3.9  ±  0.4  sec.  Thus,  for  Eurasia,  t*(f)  for  tectonic  paths  is  approximately  parallel 
to.  but  approximately  0.25  sec  (for  tp)  greater  than  t*(f)  under  the  shield.  <  ' A  /' a-  ^ 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

DO  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  DTlC  USERS  UNCLASSIFIED  _ _ 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

Dr.  Dean  Clauter _ (305)  494-52  _1 _ - 


DD  FORM  1473,84  MAR  B3  APR  edition  may  be  used  until  exhausted  _ SECURITY  CLASSIFICATION  OF  THIS  PAGE 


t  *  sr- 


All  Othe'  editions  are  obsolete 

,  -f  OZ  " 


UNCLASSIFIED 

f 


UNCLASSIFIED 


For  North  America,  the  amplitudes  of  long  period  multiple  Sc S  phases  were  analyzed  for  two 
type  of  paths.  Multiple  ScS  amplitudes  recorded  in  the  Northeastern  US  for  events  in  the  Aleutians 
(paths  across  the  Canadian  shield)  give  t&s  ~  1.1  ±  0.4  sec.  Similar  data  recorded  in  the  Southeastern 
US  for  a  Central  American  event,  gives  t scS  -  4.0  ±  0.9  sec.  •  These  data  analyses  show  that  the 
attenuation  in  the  long  period  band  in  North  America  varies  in  the  same  sense  regionally  as  shown  in 
previous  work  in  the  short  period  banu.  For  Eastern  Pacific  the  upper  mantle  attenuation  appears  to  be 
the  highest  of  all  regions  studied  with  t*  uniformly  higher  than  in  the  two  "tectonic"  continental  areas 
discussed  above. 

These  findings  support  the  general  jdea,  proven  for  continental  North  America  by  Der  and  Lees 
( 1 985)  »that  t  ,(f)  curves  for  teleseismic  paths  crossing  the  upper  mantle  in  various  areas  of  the  world 
generally  do  not  converge  towards  the  long  period  end  of  the  seismic  spectrum  as  implied  by  some 
absorption  band  shift"  models.  Instead,  the  t*(f)  curves  tend  to  run  quasi -parallel  throughout  the 
seismic  frequency  band  as  one  would  exper '  if  the  lateral  Q  variations  in  the  upper  mantle  are  caused 
by  changes  in  temperature  and  chemical  c  (position  over  large  depth  ranges  (Minster  and  Anderson, 
1980;  Jordan  1975).  These  broad-band  frequency  dependent  attenuation  models  must  be  known  when 
applying  attenuation  correction  to  yield  estimates. 

\ 

\ 


St CURITY  CLASS If ICA T IQ NjO F_T Hl$  PAGt 

UNCLASSIFIED  " 


Final  Report 


Frequency  Dependence  of  Q 


TGAL-86-5 


SUMMARY 

The  results  of  wide-band  data  analyses  for  the  determination  of  frequency  dependent  Q 
models  of  the  mantle  under  the  tectonic  areas  of  Eurasia,  North  America  and  the  Eastern 
Pacific  are  presented.  For  tectonic  areas  of  Eurasia,  the  spectra  of  short  period  P  waves  from 
Kazakh  explosions  are  used  to  estimate  t*  in  the  1  to  8  Hz  range.  For  shield  to  tectonic  paths, 
tp  -  0.35  sec  which  gives  tp  -  0.6  sec  for  tectonic  to  tectonic  paths.  P  wave  rise  times  from 
the  same  events  give  upper  bounds  on  the  estimates  of  tp  to  be  -  0.8  sec.  The  large  difference 
between  the  spectral  and  rise  time  estimates  may  be  due  to  scattering.  In  the  0.5  to  1.3  Hz 
range,  PP  and  P  amplitude  ratios  give  tp  ~  0.5  sec  for  PP’s  relative  to  surface  reflection  points 
under  the  Eurasian  shield.  Amplitude  ratios  of  long  period  multiple  ScS  phases  were  used  to 
estimate  attenuation  in  the  0.3  to  0.9  Hz  band.  Under  tectonic  Eurasia,  these  analyses  gave 
tscS  ~  3.9  ±  0.4  sec.  Thus,  for  Eurasia,  t*(f)  for  tectonic  paths  is  approximately  parallel  to,  but 
approximately  0.25  sec  (for  tp)  greater  than  t*(f)  under  the  shield. 

For  North  America,  the  amplitudes  of  long  period  multiple  ScS  phases  were  analyzed  for 
two  type  of  paths.  Multiple  ScS  amplitudes  recorded  in  the  Northeastern  US  for  events  in  the 
Aleutians  (paths  across  the  Canadian  shield)  give  tgcg  ~  1.1  ±  0.4  sec.  Similar  data  recorded 
in  the  Southeastern  US  for  a  Central  American  event  gives  tgcg  —  4.0  i  0.9  sec.  These  data 
analyses  show  that  the  attenuation  in  the  long  period  band  in  North  America  varies  in  the  same 
sense  regionally  as  shown  in  previous  work  in  the  short  period  band.  For  Eastern  Pacific  the 
upper  mantle  attenuation  appears  to  be  the  highest  of  all  regions  studied  with  t*  uniformly 
higher  than  in  the  two  "tectonic"  continental  areas  discussed  above. 

These  findings  support  the  general  idea  proven  for  continental  North  America  by  Der  and 
Lees  (1985)  that  t  (f)  curves  for  telcseismic  paths  crossing  the  upper  mantle  in  various  areas  of 
the  world  generally  do  not  converge  towards  the  long  period  end  of  the  seismic  spectrum  as 
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implied  by  some  "absorption  band  shift"  models.  Instead,  the  t*(f)  curves  tend  to  run  quasi¬ 
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in  the  upper  mantle  are  caused  by  changes  in  temperature  and  chemical  composition  over  large 
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INTRODUCTION 

This  final  report  describes  the  functional  forms  of  frequency  dependence  of  Q  ir  the  man¬ 
tle  under  various  regions  of  the  Earth.  The  study  is  exploratory  in  nature;  resolution  of  finer 
details  in  the  Q  structures  within  each  of  the  major  regions  studied  are  beyond  the  scope  of 
this  limited  project.  The  motivation  for  studies  of  this  nature  is  given  by  the  fact  that  the  form 
of  t  (f)  for  various  paths  determines  how  mb  depends  on  the  measured  values  of  the  apparent 

*  “T 

t  (t  )  derived  from  spectral  analyses.  For  some,  albeit  hypothetical,  forms  of  frequency 
dependence  of  regional  t  (f)  models,  the  relationships  between  the  measured  tp  and  mb  may  be 
quite  complex,  thus  affecting  our  ability  to  measure  attenuation-related  source  and  station 
biases  from  spectral  measurements  (Hadley  and  Mellman,  1983).  It  is  therefore  important  to 
prove  that  the  corresponding  forms  of  frequency  dependence  of  t  are  inconsistent  with  the 

data  if  we  are  to  utilize  spectral  methods  to  estimate  magnitude  bias.  Moreover,  if  the  apparent 

*  . 

t  is  not  strongly  frequency  dependent  even  though  the  absolute  t  is,  then  it  is  possible  to  esti¬ 
mate  source  region  biases  from  some  quite  simple  spectral  measurements  (De*-  et  al. ,  1979, 
1982a, b;  Der  and  Lees,  1985)  for  the  purpose  of  yield  estimation. 

The  appropriate  models  for  t  (f)  can  be  delineated  by  making  measurements  of  attenua¬ 
tion  in  different  frequency  bands  using  various  types  of  data  and  piecing  together  the  most 
likely  functional  forms  of  t  (f)  and  outlining  the  most  likely  models  of  Q(f,r)  of  the  given 
regions.  Our  measurements  in  the  long  period  band  also  have  relevance  to  acceptable  forms  of 
Q(f,r)  to  be  used  for  correcting  surface  v'1  data  for  attenuation  (Stevens,  1986;  Der,  1986). 

Various  types  of  t  (0  models  have  been  proposed  for  use  in  correcting  yield  estimates 
for  attenuation.  In  Figure  1,  we  summarize  a  few  of  these  models  for  central  and  southwestern 
United  States.  These  models  include  the  strongly  convergent  models  proposed  by  Lay  and 
Helmberger  (1981)  and  Butler  (1984),  the  weakly  convergent  models,  (that  is  less  frequency 
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dependent  than  the  earlier  pairs)  as  proposed  by  Hadley  and  Mellman  (1983)  and  the  quasi¬ 
parallel  models  by  Der  and  Less  (1986).  These  models  have  been  constructed  using  mainly 
short  period  data;  die  long  period  extrapolation  has  a  very  strong  effect  on  how  these  curves 
may  diverge  or  converge  and  on  mb  also.  It  is  the  objective  of  this  study  to  settle  such  ques¬ 
tions  for  all  the  areas  studied. 

In  studies  of  this  kind  we  are  not  looking  for  subtle  differences  in  attenuation,  but  rather 
for  changes  in  upper  mantle  Q  by  factors  exceeding  two  in  depth  ranges  exceeding  100  km.  It 
is  doubtful  that  any  lateral  changes  less  pronounced  than  this  could  be  detected  at  all  by  any 
ct  rrently  used  methodology,  especially  in  the  long  period  band  where  the  wave  amplitudes  are 
less  sensitive  to  Q.  Any  measurements  of  Q  must  be,  by  the  nature  of  the  problem,  much 
cruder  than  measurements  of  the  seismic  velocities  in  the  Earth,  for  example.  Since  the  detect¬ 
able  changes  in  t*  at  teleseismic  distances  must  be  at  least  of  the  order  of  two  and  are  mainly 
due  to  lateral  changes  of  Q  in  a  relatively  thin  upper  mantle  layer  (250-300  km),  the  Q  in  this 
layer  must  change  by  factors  considerably  larger  than  two.  Our  survey  of  the  seismological 
literature  indicates  that  some  of  the  best  documented  (indeed  classical)  studies  of  the  upper 
mantle  support  lateral  variations  in  Q  by  about  two  orders  of  magnitude  in  the  depth  range  of 
100-300  km.  Since  our  broad  band  study  encompasses  both  the  short  and  long  period  bands 
our  results  are  relevant  to  yield  estimation  from  both  short  period  body  waves  and  long  period 
surface  waves. 

In  studying  "tectonic",  or  any  regionalized  data,  a  word  of  caution  is  in  order.  Clearly, 
any  such  concept  as  a  "tectonic"  region  is  an  abstraction  that  will  have  to  be  discarded  in  the 
near  future  as  detailed  three-dimensional  velocity  models  of  the  Earth  from  tomographic  studies 
(Dziewonski  and  Anderson,  1981)  assume  concrete  forms.  In  any  case,  "shields"  correspond 
to  the  upper  mantle  regions  of  high  velocity  with  deep  "roots"  in  the  tomographic  results. 
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Tliere  is  evidence  indicating  that  tectonic  areas  are  more  heterogeneous  than  shields  with 
respect  to  Q. 


I 

R 


In  studies  of  this  sort,  the  most  difficult  task  is  to  extract  relevant  data  from  a  mass  of 
recordings  and  to  overcome  the  severe  limitations  imposed  on  the  scientist  by  geography  and 
availability  of  suitable  stations  and  events.  Expectations  with  regard  to  the  availability  of  cer¬ 
tain  kinds  of  arrivals  are  often  not  met.  In  continental  areas,  our  experience  parallels  that  of 
Sipkin  and  Jordan  (1980)  who  have  found  that  ScS  phases  are  often  masked  by  a  high  back¬ 
ground  consisting  of  continental  higher  and  fundamental  mode  waves  excited  by  the  event.  In 
our  experience,  only  a  fraction  (less  than  one  third)  of  the  events  examined  have  yielded  useful 
ScS  data  even  when  the  data  set  was  restricted  to  very  large  and  deep  events.  In  general,  one 
has  to  search  through  all  plausible  events  in  a  region  recorded  over  many  years  to  gather 
sufficient  data  for  a  study.  The  film  library  of  the  Center  for  Seismic  Studies  covers  the  years 
1970-1982  fairly  well  by  acting  as  the  main  source  of  analog  data.  Outside  this  time  period 
the  coverage  is  sporadic  and  it  was  necessary  to  acquire  data  from  other  sources  during  this 
project  (more  complete  sets  of  film  chips  exist  at  Lamont-Doherty  Geological  Observatory). 

In  previous  work  we  have  investigated  the  shield  areas  of  Eurasia  using  broad-band  t* 
measurements  and  determined  a  frequency  dependent  Q  model  for  this  region  (Der  et  al. , 
1985b).  This  model  gives  us  a  useful  basis  of  reference  for  comparing  the  frequency  depen¬ 
dence  of  Q  in  various  areas  of  the  world  chosen  for  this  project.  Following  the  methodology 
of  our  shield  studies,  this  study  avoids  mixing  global  and  regional  data  in  order  to  arrive  at  Q 
models  applicable  to  each  region. 

The  general  plan  of  this  study  is  to  make  t*  estimates  in  several  overlapping  limited  fre¬ 
quency  bands  and  combine  them  to  construct  the  best  fitting  Q(r,f)  models  that  are  the  most 
compatible  with  the  various  kinds  of  measurements.  Since  thus  far  no  convincing  evidence  has 
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surfaced  favoring  rapid  variation  in  the  derivative  of  Q  (or  t  )  with  respect  to  frequency,  the 
use  of  quasi-constant  apparent  Q  in  narrow  frequency  bands  is  not  likely  to  introduce  appreci¬ 
able  errors  in  the  construction  of  a  broad-band  Q(f)  model.  Moreover,  since  we  have  made  a 
thorough  study  of  some  shield  regions,  we  may  now  make  relative  comparisons  of  attenuation 
estimates  between  the  tectonic  regions  and  the  known  shield  regions.  There  is  also  hardly  any 
need  to  analyze  additional  short  period  data  for  North  America  in  constructing  frequency 
dependent  Q  models  for  the  region  since  it  has  been  well  studied  in  this  frequency  band  (Lay 
and  Helmberger,  1981;  Der  et  al. ,  1982).  We  will  concentrate  in  obtaining  new  long  period 
estimates  of  Q  to  supplement  the  available  short  period  estimates  from  previous  studies  to  con¬ 
struct  a  broad-band  Q(f)  model  for  North  America. 

Both  amplitude  and  spectral  measurements  are  important  in  this  study  of  Q.  In  our  data 

analyses  we  make  a  distinction  between  results  for  the  absolute  t*  and  the  apparent  t*  defined 
+  * 

as  t  .  At  a  constant  frequency,  t  can  be  determined  from  the  change  in  amplitude  spectrum  of 
an  arrival, 

A  =  e-*"',  where  t‘  =  J ^travel  time)  (1) 

Waveforms  and  spectral  ratios  are  mostly  shaped  by  the  apparent  t*  defined  as  t\  which  can  be 
written  as: 
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1980,  1982a).  This  relationship  can  best  be  tested  in  the  0.1  to  2  Hz  frequency  range  where 
attenuative  effects  are  observable  in  both  P  and  S  waves. 

In  studies  of  this  kind  the  most  difficult  part  of  the  work  is  the  finding  of  suitable  data. 
When  searching  for  data,  certain  paths  are  hard  to  access  due  to  the  location  of  the  available 
earthquakes  and  geophysical  and  political  constraints.  We  will  be  investigating  all  possible  data 
analysis  schemes  to  process  the  available  data  to  obtain  attenuation  measurements  covering  as 
broad  a  frequency  band  as  possible. 

In  the  first  part  of  this  report  we  shall  describe  the  results  obtained  for  multiple  ScS 
waves.  This  is  followed  by  our  data  analyses  for  all  of  our  individual  regions  and  some  sur¬ 
rounding  but  relevant  areas.  The  report  is  concluded  by  defining  the  most  likely  forms  of  t  (f) 
for  each  individual  region  and  present  our  conclusions  as  to  what  the  results  imply  for  yield 
estimation. 

The  reason  for  this  manner  of  presentation  is  that  the  methodology  for  determining  Q 
from  multiple  ScS  data  by  various  workers  is  quite  varied  and  the  results  are  not  directly  com¬ 
parable,  and  we  thus  have  decided  that  the  regional  variations  in  the  QScs  should  be  presented 
first  using  a  uniform  methodology  for  its  estimation.  By  bringing  all  of  our  ScS  results 
together,  as  well  as  reworking  some  of  the  data  from  other  authors,  we  are  able  to  define 
regional  patterns  and  compare  our  results  to  other  work  within  the  framework  of  a  common 
methodology.  We  follow  this  by  an  extension  of  the  results  to  shorter  periods  by  combining  a 
variety  of  new  data  with  known  regional  attenuation  patterns  for  most  of  the  regions  included 
in  this  study  (Der  ct  al. ,  1982b;  Lay  and  Helmberger,  1981).  All  data  spanning  both  long  and 
short  period  bands  are  then  combined  to  construct  the  most  likely  t*(f)  curves  for  the  regions 
of  interest. 
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ANALYSES  OF  MULTIPLE  ScS  DATA  IN  ALL  REGIONS 


Background 

The  quality  factors  derived  from  the  analysis  of  multiple  ScS  provide  a  sampling  of  the 
whole  mantle.  Since,  presumably,  most  of  the  attenuation  probably  occurs  in  the  upper  man¬ 
tle,  the  measurement  of  QScS  is  likely  to  be  the  most  indicative  of  the  dissipation  of  energy 
within  the  upper  mantle  low  Q  zone.  Although  some  researchers  postulated  a  low  Q  region 
near  the  core-mantle  boundary  (Anderson  and  Given,  1982),  the  existence  of  such  feature  is 
still  being  debated.  Prominent  multiple  ScS  phases  with  impulsive  signatures  are  frequently 
present  on  recordings  of  earthquakes  with  hypocenters  below  the  lithosphere,  these  arrivals  are 
obscured  by  high  amplitude  surface  waves  for  shallower  events.  The  SH  component  of  these 
arrivals  is  eminently  suited  for  attenuation  studies  because  of  the  lack  of  S-P  conversions  at  the 
free  surface  or  the  core  mantle  boundary.  This  study  analyzes  multiple  ScS  phases  in  the  time 
domain  to  obtain  QScS  for  different  paths  across  the  tectonic  S.  Pacific,  continental  S.  America 
and  N.  America  and,  in  doing  this,  to  detect  any  regional  variations  of  QScS. 

Various  authors  have  utilized  ScS  arrivals  to  obtain  average  Q  values  for  the  whole  man¬ 
tle.  The  methodologies  applied  vary  considerably.  Jordan  and  Sipkin  (1977)  applied  a  stack¬ 
ing  method  in  the  spectral  domain  to  obtain  a  statistical  measurement  of  QScS  for  different 
parts  of  the  world.  Nakanishi  (1979)  obtained  QScS  for  the  Sea  of  Okhotsk,  Sea  of  Japan  and 
the  lzu-Bonin  region  using  a  maximum  likelihood  estimates  on  spectral  ratios  between  multiple 
ScS  phases  from  seismic  recordings  in  Japan.  Sato  and  Espinosa  (1967)  obtained  a  lower 
bound  lor  QScS  for  paths  from  S.  America  to  N.  American  stations  using  essentially  a  spectral 
ratio  technique.  They  also  obtained  estimates  of  the  core  viscosity  and  postulated  the  frequency 
dependence  of  Q.  Similar  finding  was  reported  by  Yoshida  and  Tsujiura  (1975).  It  is  quite 
difficult  to  compare  the  results  obtained  by  the  various  authors  since  all  methods  have  various 
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degrees  of  inherent  bias.  The  method  used  to  estimate  Q  from  multiple  ScS  is  largely  a  matter 
of  personal  preference  and  philosophy. 

Estimates  of  QScS  from  the  shapes  of  the  spectra  as  opposed  to  the  decay  of  absolute 
amplitudes  in  the  successive  ScS  arrivals  are  heavily  influenced  by  the  ever-present  scattered 
surface  wave  background  on  the  long  period  seismograms.  The  ScS  pulses  are  of  very  short 
duration  thus  making  the  window-carpentry  for  spectral  analyses  very  difficult.  Such  noise 
contamination  will  result  in  overestimates  of  Q.  Spectral  stacking  methods  may  de-emphasize 
the  high  frequency  content  of  the  data  thus  biasing  the  Q  values  to  be  lower. 

We  have  used  two  time  domain  approaches  to  estimate  Qscs  for  different  regions  of  the 
world.  A  robust  measurement  of  Qscs  is  introduced  in  this  study  using  a  time  domain 
waveform  matching  scheme  to  obtain  the  Q  operator  that  produces  both  the  broadening  and 
decrease  in  amplitude  of  the  successive  multiply  reflected  ScS  phases.  By  processing  a  large 
amount  of  data  from  various  parts  of  the  world  we  obtained  Q  estimates  that  are  comparable 
because  of  the  identical  processing  method  applied.  By  reprocessing  the  data  from  some  ear¬ 
lier  analyses  we  are  also  able  to  compare  various  regions  of  the  earth  in  a  common  framework. 

The  other  approach  has  been  used  in  our  earlier  work  on  Q  under  the  Eurasian  shield 
(Der  et  al.,  1984)  utilizing  WWSSN  photographic  records  to  obtain  QScS.  For  these  we  used 
a  simplified  method  of  analysis  by  fitting  the  decay  rate  of  the  successive  ScS  pulses  (corrected 
for  geometrical  spreading  effects)  graphically  to  values  of  t  .  In  such  analyses  we  assumed 
that  the  ScS  pulses  have  a  dominant  frequency  and  that  the  pulses  are  monochromatic.  This 
method  was  found  to  be  generally  compatible  with  the  more  precise  time  domain  fitting 
methods  applied  to  our  digital  data.  This  method  is  mainly  used  for  N.  America  and  Eurasia 
where  the  station  and  earthquake  configurations  prevent  us  to  obtain  adequate  digital  data  to 
perform  the  time  domain  waveform  matching  scheme.  Results  on  t  lor  North  America  and 


I’eledyne  Geo  tec  It 


8 


October  1986 


Final  Report 


Frequency  Dependence  of  Q 


TGAL-86-5 


Eurasia  are  mostly  obtained  using  this  method.  Although  hand-digitization  was  attempted  on 
some  of  die  photographically  recorded  data,  the  lack  of  precision  in  lining  up  the  horizontal 

components  for  rotation,  prevented  us  to  gain  any  more  precision  compared  to  the  simplified 
analysis. 

Methodology 

The  time  domain  waveform  matching  scheme  compares  the  entire  waveform  of  ScS  with 
all  of  its  multiples.  We  convolve  the  first  well  recorded,  clear  ScSj-j  phase  with  succession  of 
t  operators  such  that  the  convolved  phase  will  have  the  same  amplitude  decay  rate  and 
waveform  characteristics  as  those  of  the  succeeding  multiples.  These  ScSN  phases  are 
identified  with  the  aid  of  the  J-B  travel  time  table.  After  correcting  for  geometrical  spreading 
effects,  we  add  for  each  increase  of  multiplicity  a  trial  value  for  the  contribution  to  the  t+  for  a 
single  passage  through  the  mantle  multiplied  by  a  factor  of  two.  The  observed  amplitude  of 
each  convolved  phase  is  compared  to  the  observation  and  the  trial  value  is  changed  until  a 
good  fit  is  obtained.  Due  to  the  low  amplitude  of  the  last  observable  ScS  multiples,  comparis¬ 
ons  performed  to  these  phases  act  mainly  as  constraints  on  the  decay  rate.  Most  of  the  ampli¬ 
tude  and  phase  information  is  utilized  in  the  analysis  if  the  waveforms  fit  also,  which  is  usually 
the  case.  Examples  for  two  circum-Pacific  deep  focus  events  are  shown  in  Figure  2  to  indicate 
die  ease  of  the  method. 

This  method  has  the  advantage  over  spectral  shape  measurements  in  being  simple  and 
yields  much  more  stable  measurements.  We  feel  that  the  spectral  shapes  and  the  waveforms 
are  heavily  influenced  by  the  scattered  background  and  any  Q  values  that  conflict  with  the 
gross  amplitudes  of  successive  ScS  arrivals  should  be  considered  suspicious.  Often  we  are 
able  to  put  upper  limits  to  Q  by  measuring  the  background  level  at  the  expected  arrival  time  of 
the  next  Sc Sn  pulse  at  ter  the  last  observable  arrival.  ScS  pulses  and  their  associated  surface 
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reflections  often  form  characteristic  pairs  which  can  be  recognized  even  at  low  amplitudes. 
Some  of  the  examples  of  such  ScSn-sScSn  pairs  are  shown  in  Figure  3. 

Data  Analysis 

We  have  collected  data  from  a  variety  of  tectonic  environments  around  the  world.  Multi¬ 
ple  ScS  phases  are  the  best  seen  for  deep  earthquakes  for  which  the  surface  waves  that  com¬ 
monly  mask  ScS  phases  are  not  strongly  excited.  Event  parameters  for  all  the  earthquakes  used 
are  given  in  Table  1.  The  first  region  we  shall  study  is  the  E.  Pacific  Ocean.  The  tectonic 
subdivision  of  the  Pacific  Ocean  is  shown  in  Figure  4.  Seven  earthquakes  are  used  for  the  E. 
Pacific  path  (Figure  5).  They  are  mostly  deep  focus  earthquakes  associated  with  the  Tonga-Fiji 
trench  and  recorded  at  digital  stations  ZOBO  and  BOCO  in  S.  America.  The  bounce  points  of 
these  events  are  also  shown  in  Figure  5.  From  the  figure,  it  is  observed  that  the  bounce  points 
are  located  away  from  the  E.  Pacific  Rise  where  active  spreading  and  upwelling  of  magma  are 
occurring.  Any  average  of  t*  including  measurement  of  the  upper  mantle  beneath  or  very  near 
E.  Pacific  Rise  may  be  biased  towards  low  Q.  The  phases  used  for  these  events  are  mostly 
ScS2,  ScSy,  SCS4.  The  first  ScS  is  not  used  in  this  part  of  the  study  since  it  is  either  not 
observed  (shadow  zone)  or  arrives  too  close  to  the  S  phase  in  time.  The  waveforms  and  th  ir 
matchings  are  shown  in  Appendix  A.  In  most  cases,  the  amplitudes  of  the  multiple  ScS  decay 
quite  rapidly  as  the  number  of  bounces  increases  as  shown  in  Figure  3  and  other  examples  in 
Appendix  A.  The  values  of  tscS  across  the  S.  Pacific  are  generally  high  and  average  about  7.0 
sec. 

To  constrain  our  measurements,  we  have  also  applied  the  same  method  to  other  regions 
and  some  published  data.  Eight  events  in  the  Tonga-Fiji  trench  recorded  at  digital  stations, 
A  FI,  HON,  CTA  and  GUMO  (Figure  6,  Table  1)  have  been  used  to  study  the  QScS  for  the 
tectonic  S.W.  Pacific.  These  phases  are  ScS,  ScSn,  ScSy  and  ScS4  with  surface  reflection 
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ri^ui-c  3.  Examples  of  seismograms  showing  the  decay  rate  of  ScSN  with  multiples.  ScS 

multiples  lire  observable  up  to  ScS4  for  shield  S.  America  and  W.  Pacific  paths  whereas 
ScS4  is  into  the  noise  for  path  crossing  the  E.  Pacific. 
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Table  1 


Parameters  of  earthquakes  used  in  ScS  study  of  tectonic  Pacific 


Path:  Fiji-Tonga  events  recorded  at  ZOBO  and  BOCO 


Event  # 

date 

latitude 

longitude 

depth 

mb 

1 

Nov  25,  1976 

-19.50 

-177.70 

424 

2 

Jul  17,  1978 

-23.63 

-179.96 

519 

5.4 

3 

Apr  24,  1979 

-20.82 

-178.67 

596 

5.9 

4 

Jul  20,  1980 

-17.88 

-178.61 

596 

5.9 

5 

Sep  28,  1981 

-29.33 

-179.04 

323 

6 

Nov  25,  1981 

-20.94 

-178.86 

614 

5.9 

7 

Sep  16,  1983 

-24.03 

-179.80 

Path:  Fiji-Tonga  events  recorded  at  AFI 

,  HON  and  GUMO 

Event  # 

date 

latitude 

longitude 

depth 

mb 

8 

Sep  28,  1981 

-29.33 

-179.04 

323 

9 

Jan  19,  1984 

-23.64 

-178.32 

332 

5.8 

10 

May  30,  1984 

-04.84 

151.58 

174 

6.2 

11 

Aug  26,  1984 

-23.59 

179.07 

560 

5.9 

12 

Nov  17,  1984 

-18.79 

-178.03 

451 

6.1 

13 

Nov  22,  1984 

-17.78 

-178.05 

646 

5.9 

14 

Dec  12,  1985 

-21.66 

-176.38 

155 

5.9 

15 

Aug  28,  1985 

-21.01 

-178.98 

625 

6.1 

Path:  South  American  events  recorded  at  ZOBO  and  BDF 

Event  # 

date 

latitude 

longitude 

depth 

mb 

16 

Feb  04,  1977 

-22.66 

-63.39 

555 

5.9 

17 

May  21,  1979 

-15.44 

-70.05 

209 

5.9 

18 

Jan  10,  1983 

-22.01 

-68.47 

121 

5.6 

19 

Jun  02,  1983 

-09.51 

-71.25 

599 

5.9 

20 

Dec  21,  1983 

-28.19 

-63.17 

602 

6.2 

21 

May  01,  1985 

-09.20 

-71.23 

600 

22 

Jun  10,  1985 

-27.96 

-66.99 

151 

5.8 

23 

Nov  29,  1985 

-22.76 

-63.62 

521 

5.3 

Path:  Events  used  by  Nakanishi  (1979)  for  NW  Pacific 

Event  # 

date 

latitude 

longitude 

depth 

mb 

24 

Jun  26,  1975 

+31.90 

+  138.20 

394 

5.3 

25 

Dec  21,  1975 

+51.90 

+  151.60 

546 

26 

Aug  12,  1975 

+32.10 

+  137.90 

399 

5.6 

27 

Jul  28,  t973 

+50.50 

+  148.90 

585 

5.5 

Path:  Events  used  by  Sipkin  and  Jordan 

(1977)  for  die  Pacific 

Event  # 

date 

latitude 

longitude 

depth 

mb 

28 

Apr  10,  1976 

-17.10 

-178.46 

548 

5.7 

29 

Nov  25,  1976 

-19.50 

-177.70 

424 

30 

Jan  9,  1976 

-15.80 

+  167.85 

173 

5.7 

K 
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points  within  the  S.W.  Pacific.  The  seismograms  and  waveform  matches  are  presented  in 
Appendix  A.  The  resulting  average  tscs  for  this  relatively  limited  region  is  about  4.0  sec. 

We  have  also  applied  the  method  to  shield  S.  America.  Deep  focus  S.  American  trench 

earthquakes  recorded  at  digital  stations  ZOBO  and  BDF  with  bounce  points  within  the  region 

are  shown  in  Figure  7  and  listed  in  Table  1.  The  waveforms  and  data  analysis  runs  are  shown 

in  Appendix  A.  The  tscs  for  the  shield  S.  Ameiica  is  around  2.0  seconds  whereas  for  the  tec- 
•  * 

tonic  paths,  tscs  is  up  to  5.0  sec. 

We  have  also  re-analyzed  some  data  used  by  Nakanishi(1979)  and  Sipkin  and  Jor- 
dan(1980).  These  include  Izu-Bonin  and  Kurile  events  recorded  at  HGLP  digital  station  MAT 
with  ScS,  ScS2,  SCS3  and  SCS4  sampling  the  NW  Pacific.  Figure  8  shows  the  locations  of 
stations  and  reflection  points  with  the  seismograms  shown  in  Appendix  A.  The  tscS  obtained 
for  this  region  is  around  4.0  sec  and  being  very  close  to  Nakanishi  (1979)  results.  Several 
paths  studied  by  Sipkin  and  Jordan  (1980)  are  reprocessed  here  and  the  locations  of  the  surface 
reflection  points  are  shown  in  Figure  9.  Two  paths  are  shown  in  Figure  9a  sampling  the  W. 
Pacific  and  one  path  with  ScS^  bounces  across  the  Pacific  are  shown  in  Figure  9b.  The 
waveforms  and  data  analysis  runs  are  shown  in  Appendix  A.  Most  results  are  similar,  but  the 
values  derived  by  Sipkin  and  Jordan  (1980)  tend  to  be  somewhat  higher  than  ours.  All  ts*cs 
results  are  summarized  in  Table  2  and  are  compared  to  those  obtained  by  the  amplitude  decay 
method  for  Eurasia  and  N.  America  as  summarized  below. 

Multiple  ScS  data  were  also  collected  for  deep  Hindu  Kush  earthquakes  observed  at  vari¬ 
ous  WWSSN  station.  Figure  10  shows  the  tectonic  and  shield  regions  of  Eurasia.  These  paths 
sample  two  of  our  chosen  regions  fairly  well  as  indicated  by  the  map  in  Figure  1  1.  The  event 
parameters  are  listed  in  Tables  3  and  4.  We  have  presented  some  examples  of  the  ScS  and 
sScS  multiples  in  Figures  12-14.  The  plots  of  ScS  phases,  normalized  to  ScS2  and  plotted 
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Figure  7.  Map  showing  the  location  of  events  and  stations  in  S.  America  and  the  associated 
bounce  points  for  all  the  ScSN  phases  used  are  shown  in  +’s  for  ScS2  phases,  and  X’s  fc 
ScSv 


VV 


Wo  PACIFIC 


Figure  9a.  Map  showing  the  location  of  events  and  stations  in  the  Western  Pacific  as  used 
by  Sipkin  and  Jordan  (1977).  The  bounce  points  for  all  the  ScSN  phases  used  are  shown  as 
+  ’s  for  ScS2,  X’s  for  ScS3  and  Z’s  for  Sc S4. 
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Table  2 

A  summary  of  ScSN  phases  used  and  ts*cS  for  all  regions 


Region:  E.  Pacific 


Event  tf 


Station 


ZOBO 

ZOBO 

BOCO 

BOCO 

BOCO 

BOCO 

BOCO 

ZOBO 


Observed  phase  type 
ScS2  ScS3 


Region:  SW  Pacific 


Event  ft  Station 


A  FI 
A  FI 
HON 
GUMO 
AFI 
AFI 
HON 
AFI 
HON 
CTA 


Observed  phase  type 
ScS2  ScS3 


4.0 

5.0 

5.0 

5.0 

5.5 

2.5 
7.0 
5.0 
2.5 

5.0  ft  (6.6) 


Region:  Shield  S.  America 


Event  ft 


Station 


ZOBO 

ZOBO 

ZOBO 

ZOBO 

BOCO 

BDF 

BDF 

ZOBO 

ZOBO 


Observed  phase  type 
Sc S2  ScS3 


of  -v*  ’wrw  -j  rj  -..« . 


Region:  NVV  Pacific 


Table  2  (cont’d) 


Event  It  Station 


Observed  phase  type 
ScS2  SCS3 


4.5  *  (1.5) 
4.0  *  (3.8) 
4.0  *  (3.4) 
5.0  *  (3.5) 

2.5  #  (5.8> 


Region:  N.  Trans-Pacific 


Event  #  Station 


Observed  phase  type 


2.5  #  (5.8) 


Region:  Shield  Eurasia 


Region:  Tectonic  Eurasia 


Region:  Shield  North  America 


Region:  Tectonic  N.  America 


II  (  )  value  front  Sipkin  and  Jordan  (1977) 


3.0-3.9 


1. 1-2.5 


4  0  4.5 


*  (  )  value  front  Nakanishi  (1979) 


Table  3 


Event  Parameters  for  Far  Eastern  Events 
Used  in  Estimating  t§cS  Under  the  Eurasian  Shield 


date 

origin  time 

latitude 

longitude 

depth 

mb 

(sec) 

(degrees) 

(degrees) 

(km) 

(ISC) 

660606 

07:46:16.1 

36.4 

71.1 

221.0 

6.2 

731017 

03:16:18.6 

36.4 

71.2 

221.0 

5.5 

740730 

05:12:40.6 

36.4 

70.8 

211.0 

6.5 

780421 

15:22:57.9 

36.63 

71.26 

229.6 

5.8 

790816 

21:31:26.3 

41.81 

130.19 

588.0 

6.1 

810502 

16:05:02.2 

36.8 

70.99 

229.0 

5.9 

s 


Table  4 


Event  Parameters  for  Hindu  Kush  Events 
Used  in  Estimating  t*ScS  Under  Eurasian  Tectonic  Areas 


date 

origin  time 

latitude 

longitude 

depth 

mb 

(sec) 

(degrees) 

(degrees) 

(km) 

(ISC) 

740730 

05:12:40.6 

36.35 

70.76 

211.0 

6.5 

761127 

21:42:12.1 

36.52 

71.04 

190.2 

6.2 

790820 

03:50:52.3 

36.48 

70.15 

230.5 

5.9 

60  sec 


icings  of  multiple  ScS  phases  followed  by  the  corresponding  sScS  multiples 
;tation  JER  for  the  Hindu-Kush  earthquake  of  8/20/79. 
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against  f(N-2)  where  f  is  frequency  and  N  is  the  order  of  ScS  phase  are  shown  in  Figure  15 
for  the  two  source  regions.  The  resulting  fit  to  the  data  gives  a  t<;c«;  of  2.5  sec  for  shield 
Eurasia  and  3.9  sec  for  tectonic  Eurasia. 

For  North  America,  the  nearby  deep  earthquake  sources  are  Central  America  and  the 
Alaskan-Aleutian  region.  Stations  with  multiple  ScS  phase  surface  reflection  point  in  the  tec¬ 
tonic  and  shield  regions  are  found  in  the  western  and  eastern  United  States  respectively.  The 
tectonic  and  shield  boundary  in  N.  America  is  shown  in  Figure  16.  We  found  it  difficult  to 
obtain  events  in  Central  America  with  multiple  ScS  phases  as  spectacular  as  those  in  Eurasia 
and  Pacific.  The  deepest  earthquakes  in  this  area  are  relatively  shallow,  not  much  more  than 
100  km  in  depth  and  thus  surface  waves  tend  to  obscure  ScS  arrivals.  Moreover,  the  event 
magnitudes  are  relatively  low  compared  to  other  regions  with  deep  earthquakes,  the  seismicity 
is  also  relatively  low. 

In  the  Alaskan-Aleutian  region  the  sources  are  somewhat  deeper  and  the  seismicity  higher 
but  events  towards  the  western  end  of  the  Aleutian  chain  often  have  surface  reflection  points 
for  multiple  ScS  beneath  the  ocean.  The  events  we  judged  to  be  suitable  for  analysis  are  listed 
in  Tables  5  and  6  and  their  locations  with  the  recording  stations  and  the  surface  reflection 
points  for  ScS  phases  are  shown  in  Figure  17. 

Figures  18-20  show  some  tracings  of  multiple  ScS  data  for  the  two  source  regions.  In 
general,  as  seen  in  these  figures,  the  ScS  phases  for  events  in  Central  America  with  bounce 
points  within  the  western  US  fall  off  with  increasing  order  much  faster  than  those  reflecting 
under  die  Canadian  shield.  Moreover,  for  some  of  the  ScS  phases  with  Canadian  paths,  the 
period  is  quite  short  and  does  not  change  much  with  increasing  multiplicity  of  ScS.  These 
facts  indicate  very  low  attenuation  under  Canada  and  much  higher  attenuation  under  the  tecton¬ 
ically  active  regions  of  Mexico  and  southwestern  United  States. 
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Table  5 


Event  Parameters  for  Aleutian  Events  Used  in  Estimating 

t? 


t*Sc s 

Under  North 

America 

date 

origin  time 

latitude 

longitude 

depth 

mb 

(sec) 

(degrees) 

(degrees) 

(km) 

(ISC) 

)0228 

10:52:31.2 

52.70 

-175.10 

162. 

6.1 

.0405 

09:04:42.8 

53.36 

-170.55 

153. 

5.8 

0904 

15:53:25.4 

54.98 

-163.36 

107. 

5.8 

m. 


Table  6 


Event  Parameters  for  Central  American  Event  Used 
in  Estimating  t*scs  Under  North  America 


date 

origin  time 

latitude 

longitude 

depth 

mb 

(sec) 

(degrees) 

(degrees) 

(km) 

(ISC) 

740306 

01:40:26.4 

12.29 

-86.39 

110. 

5.8 

3G 


720421 


WES  EW 


Sc  sScS^ 

1  1 


1 _ 1 

60  sec 


Figure  18.  Tracings  of  multiple  ScS  arrivals  with  the  corresponding  sScS  phases  at  the 
WWSSN  station  WES.  These  phases  have  relatively  short  periods  and  low  decay  rates  thus 
indicating  a  high  Q  path. 


710405 

BLA 

NS 

sS  ScS 

sScS 
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I  I  I 


ScS2  sScS2 

l  i 

-'V* 


Sc S3  SSCS3 

1  1 


60  sec 


Figure  19.  Tracings  of  multiple  ScS  arrivals  with  the  corresponding  sScS  phases  from  the 
4/05/71  Aleutian  event  as  recorded  at  the  WWSSN  station  BLA. 


740306  TUC  NS 


ScS- 


•  ^  v-;;  5? 


60  sec 


Figure  20.  Tracings  of  multiple  ScS  arrivals  with  the  corresponding  sScS  phases  for  the 
3/06/74  Central  American  earthquake  as  recorded  at  the  WWSSN  station  TUC.  These  phases 
have  a  rapidly  decreasing  amplitude  due  to  high  attenuation  in  the  mantle  under  Mexico. 
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The  plots  of  ScS  phases,  normalized  to  ScSi  and  plotted  against  f(N-2)  where  f  is  fre¬ 
quency  and  N  is  the  order  of  ScS  phase  are  shown  in  Figure  21  for  the  two  source  regions. 

These  data  give  ts*cS  values  of  1.1  sec  for  a  single  ScS  leg  under  Canada  and  4.0  sec  under 

( 

Mexico-Arizona.  It  appears  that  the  regional  Q  differentials  are  no  less  significant  and  are  of 
the  same  sign  in  the  long  period  band  as  in  the  short  period  band  in  accordance  with  the 
hypothesis  of  quasi-parallel  t*(f)  curves  (Dcr  and  Lees,  1985).  The  results  are  summarized  in 
Table  2. 

In  summary,  the  regional  differentials  in  attenuation  are  similar  in  sense  that  those 
reported  in  earlier  studies.  We  have  obtained  an  average  ts*cS  of  7.0  sec  for  paths  across  the  E. 
Pacific  while  the  SW  Pacific  and  the  NW  Pacific  are  characterized  by  moderate  attenuation 
with  t.ScS  values  close  to  4.0  sec.  The  results  for  the  Eurasian,  North  and  South  American 
shield  areas  are  significantly  lower  in  ts*cS  values,  mostly  in  the  2-3  sec  range.  The  tectonically 
active  continental  areas  of  Eurasia  and  North  America  are  characterized  by  higher  t<;cS  esti¬ 
mates  typically  near  4.0  sec. 

Discussion  of  the  ScS  Results 

The  range  of  variations  in  the  average  tscS  contributed  by  the  passage  of  SH  waves 
through  the  upper  mantle  seems  to  be  at  least  3.5:1.  Although  this  range  may  be  increased  by 
die  inevitable  scatter  in  any  t<;cS  estimation  procedures,  the  differences  are  significant.  The 
measurement  of  mantle  averaged  Q  values  does  not  provide  enough  information  about  the 
forms  of  the  depth  dependence  of  Q  in  the  regions  studied.  It  is  reasonable  to  assume,  how¬ 
ever,  that  most  of  the  lateral  variation  in  the  mantle  Q  is  associated  with  the  lateral  changes  in 
die  properties  of  the  upper  mantle  within  a  zone  200-300  km  thick,  since  this  is  the  vertical 
extent  of  the  largest  lateral  variations  of  the  seismic  velocities  as  well  (Dziewonski  and  Ander¬ 
son,  1981;  Grand  and  Heimberger,  1984;  Jordan,  1975),  presumably  due  to  a  common  cause, 
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lateral  changes  in  temperature  or  ehemical  composition.  The  large  variations  in  the  mantle 
averaged  Q  values  imply  much  larger  proportional  ehanges  in  the  upper  mantle  Q  in  this  layer. 

Our  analyses  of  multiple  SeS  data  show,  similarly  to  other  previous  studies,  that  the 
attenuation  properties  in  the  mantle  vary  considerably  from  region  to  region.  This  variation  is 
similar  to  the  corresponding  variations  in  the  short  period  band  (Der  et  al.,  1982b)  in  sign. 
Our  Q  values  vary  similarly  to  those  reported  by  Sipkin  and  Jordan  (1980)  but  are  significantly 
higher  than  theirs  as  shown  in  Table  2.  From  ScS  data  alone  we  eannot,  of  eourse,  delineate  a 
complete  profile  of  the  depth  dependence  of  Q,  nor  any  frequency  dependence  of  Q  in  the  nar¬ 
row  frequency  band  we  have  been  working  in. 

The  results  are  important  with  regards  to  possible  is*cS(f)  models.  Assuming  that  the 
lateral  variations  in  Q  we  observe  with  the  multiple  ScS  data  are  confined  in  a  200-250  km 
thiek  layer  in  the  upper  mantle,  wherever  the  lateral  velocity  variations  are  the  most  pro¬ 
nounced,  presumably  due  to  the  same  causes,  temperature  and  chemical  variations  in  the  upper 
mantle,  (Jordan,  1975)  these  results  imply  lateral  variations  in  of  up  to  an  order  of  mag¬ 
nitude. 
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ANALYSES  OF  DATA  FROM  THE  TECTONIC  AREAS  OF  EURASIA 

Analyses  of  I*  Wave  Spectra  from  Nuclear  Explosions 

Analyses  of  P  wave  spectra  of  nuclear  explosions  in  Kazakh  consist  of  estimating  the 
best  fitting  tp  values  from  the  spectra.  The  instrument  response  and  a  von  Seggern-Blandford 
(1972)  explosion  source  are  removed  from  die  spectra  of  each  P  wave  and  the  slope  of  a  linear 
le  uC  squares  fit  line  through  the  reduced  spectrum  within  the  region  of  good  signal-to-noise 
ratio  provides  an  estimate  of  tp.  Sample  P  wave  spectra  are  shown  in  Figure  22  for  the  tec¬ 
tonic  stations  MAIO  and  ANTO  and  the  shield  station  KONO.  The  average  tp  obtained  for 
the  shield-tectonic  path  is  around  0.28  sec.  The  tectonic  stations  show  a  much  more  limited 
range  of  good  signal  to  noise  than  does  KONO.  Except  for  TATO,  at  least  10  observations 
were  used  for  each  station.  The  event  parameters  for  the  explosions  used  in  this  study  are  in 
Table  7  and  the  values  obtained  are  given  in  Table  8. 

The  values  in  Table  8  are  not  dissimilar  to  those  quoted  by  Der  et  al.  (19S5b)  and  are 
not  particularly  high.  They  are  definitely  less  than  t*  values  for  paths  involving  the  upper  man¬ 
tle  under  the  Basin  and  Range  (Der  et  al.,  1985b;  Bache  et  al.,  1985).  Since  t+  values  for 
paths  between  Kazakh  and  Scandinavia  are  near  0.15  sec,  only  values  for  the  stations  ANTO, 
CHTO,  GRFO  and  MAIO  can  be  called  higher  than  the  typical  values  for  shields.  Since  the 
sources  of  these  events  were  in  high  Q  area,  the  additional  attenuation  must  be  attributed  to 
attenuation  in  the  mantle  under  these  stations. 

These  findings  are  in  agreement  with  some  results  of  one  of  our  previous  studies  (Der 
et  al.,  1982b)  in  which  we  found  that  unlike  the  shields,  the  so  called  "tectonic"  areas  of  the 
world  are  quite  heterogeneous  with  respect  to  Q  and  both  high  and  low  Q  paths  are  possible 
through  such  regions.  This  manifests  itself  in  the  large  scatter  in  the  t*  estimates  for  paths 
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Figure  22a.  Samples  of  P  waves  spectra  from  Kazakh  nuclear  explosions,  t-’s  are  the  signal 
and  the  o’s  are  the  noise  level  for  an  event  1976241  recorded  at  MAIO. 
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Figure  22b.  Samples  of  P  waves  spectra  from  Kazakh  nuclear  explosions.  +’s  are  the  signal 
and  the  o’s  are  the  noise  level  for  an  event  1979301  recorded  at  ANTO. 
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Figure  22c.  Samples  of  P  waves  spectra  from  Kazakh  nuclear  explosions.  +’s  are  the  signal 
and  the  o’s  are  the  noise  level  for  an  event  1979230  recorded  at  KONO. 


Table  7 


Event  Parameters  for  Kazakh  Explosions 


■/ 

* 


date 

origin  time 

latitude 

longitude 

depth 

mb 

(sec) 

(degrees) 

(degrees) 

(km) 

(ISC) 

1976186 

2:56:57.480 

49.85 

78.97 

0.00 

5.80 

1976241 

2:56:57.600 

49.95 

78.98 

0.00 

5.80 

1976328 

5:02:57.470 

49.97 

79.01 

0.00 

5.80 

1976342 

4:56:57.480 

49.87 

78.89 

0.00 

5.90 

1977149 

2:56:57.490 

49.85 

78.84 

0.00 

5.80 

1977248 

3:02:57.760 

50.05 

78.93 

0.00 

5.80 

1977334 

4:06:57.630 

49.93 

78.89 

0.00 

6.00 

1978085 

3:56:57.690 

49.71 

78.06 

0.00 

5.60 

1978162 

2:56:57.800 

49.88 

78.81 

0.00 

5.90 

1978186 

2:46:57.510 

49.84 

78.91 

0.00 

5.80 

1978209 

2:46:57.760 

49.73 

78.15 

0.00 

5.70 

1978241 

2:37:06.380 

49.98 

79.02 

0.00 

5.90 

1978258 

2:36:57.500 

49.91 

78.94 

0.00 

6.00 

1978308 

5:05:57.550 

50.03 

78.97 

0.00 

5.60 

1978333 

4:33:02.870 

49.93 

78.77 

0.00 

6.00 

1979174 

2:56:59.050 

49.88 

78.92 

0.00 

6.20 

1979188 

3:46:57.480 

50.05 

79.06 

0.00 

5.80 

1979216 

3:56:57.250 

49.86 

78.94 

0.00 

6.10 

1979230 

2:51:57.310 

49.93 

78.98 

0.00 

6.10 

1979301 

3:16:57.030 

49.96 

79.07 

0.00 

6.00 

1979336 

4:36:57.620 

49.88 

78.83 

0.00 

6.00 

1979357 

4:56:57.640 

49.92 

78.80 

0.00 

6.20 

1980116 

3:56:57.440 

49.92 

78.81 

0.00 

5.50 

1980143 

3:56:57.780 

49.75 

78.11 

0.00 

5.50 

1980164 

3:26:57.670 

49.95 

79.05 

0.00 

5.60 

1980181 

2:32:57.750 

49.90 

78.86 

0.00 

5.70 

1980258 

2:42:39.300 

49.98 

78.89 

0.00 

6.20 

1980286 

3:34:14.100 

49.96 

79.08 

0.00 

5.90 

1980349 

3:47:06.600 

49.93 

79.00 

0.00 

5.90 

1980362 

4:09:08.200 

50.04 

79.05 

0.00 

5.90 

1981088 

4:03:50.000 

50.01 

79.02 

0.00 

5.60 

1981112 

1:17:11.400 

49.90 

78.90 

0.00 

5.90 

■l 


Table  8 


t*|>  Estimates  for  Paths 
from  East  Kazakh  to  Eurasian  Stations 
for  Shield  -  Tectonic  Paths 


station 

frequency  range 

t* 

(Hz) 

(sec) 

ANTO 

0.5-4 

0.35±0.04 

CHTO 

0.5-5 

0.29±0.06 

GRFO 

0.5-4 

0.3410.06 

KAAO 

0.5-6 

0.1310.03 

KONO 

0.5-5 

0.2410.05 

MAIO 

0.5-4 

0.2810.04 

MAJO 

0.8-5 

0.2010.05 

SHIO 

0.5-5 

0.2210.04 

TATO 

1.0-3 

0.1710.07 
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crossing  tectonically  active  areas. 


Analyses  of  Short  Period  PP/  P  Spectral  Ratios 


Unlike  the  analyses  of  direct  P  wave  spectra  which  are  diagnostic  of  the  upper  mantle 
under  the  station  and  source,  comparisons  of  the  frequency  contents  of  P  and  PP  will  give 

information  about  the  mantle  attenuation  under  the  surface  reflection  point  of  PP. 

s 

Short  period  PP  arrivals  arc  commonly  engulfed  in  the  coda  of  the  first  P  arrival. 


Although  some  successful  analyses  of  PP/P  spectral  ratios  have  been  made  on  carefully 


selected  data  to  determine  Q  (Shore,  1983),  we  have  found  it  to  be  more  expedient  to  estimate 
Q  by  band-pass  filtering  the  wavetrain  consisting  of  P  and  PP  arrivals,  the  same  window  also 


includes  sP  for  earthquakes.  Moreover,  at  some  distances  PcP  is  often  near  these  arrivals.  It 


is  expected  that  the  PP  arrival  will  disappear  if  the  surface  reflection  point  is  in  highly  attenuat- 

i 

ing  regions,  and  the  high  frequency  expression  of  the  phase  should  disappear  first. 


We  have  processed  a  number  of  PP  and  P  pairs  with  sharply  tuned  (48  dB/octave)  band 


pass  filters  to  minimize  the  effect  of  spectral  falloffs  of  the  seismic  sources  on  the  outputs. 


The  events  used  for  the  PP  analyses  are  all  explosions  at  Kazakh.  We  chose  these  because  the 


time  functions  of  these  sources  are  relatively  short  in  duration  compared  to  earthquakes  and  PP 


is  easier  to  recognize  on  the  records.  Moreover,  explosions  generally  excite  high  frequency 


waves  more  efficiently.  The  same  approach  is  also  used  for  earthquakes  for  the  Pacific  region 


as  discussed  elsewhere  in  this  report. 


Appendix  B  shows  some  of  the  bandpass -filtered  waveforms  used  for  the  PP  study.  At 


some  stations,  such  as  KONO,  SI  110,  and  GRFO,  the  PP  phases  appear  frequently  and  can  be 


observed  continuously  as  we  go  towards  bands  with  increasing  center  frequency.  The  PcP 


phase  is  also  often  visible.  At  other  stations,  such  as  CUTO  and  ANTO,  the  PP  phase  is  either 


not  observed  at  all  or  it  disappears  into  the  noise  background  at  higher  frequencies.  The 
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locations  of  the  surface  reflection  points  together  with  the  estimated  t*  are  shown  in  Figure 
23a,  and  the  results  are  tabulated  in  Table  9.  For  some  events  at  GRFO  and  KONO,  die  PP/P 
ratio  initially  increases  with  frequency  thus  making  the  overall  t*  estimates  negative.  We  do 
not  understand  the  reasons  for  this,  but  this  is  clearly  an  indication  of  a  high  Q  mantle  under 
die  surface  reflection  point  which  lie  in  the  Russian  shield.  In  Figure  23a  we  designate  these 
with  the  t  estimate  of  zero.  Since  we  do  not  believe  in  negative  or  zero  t*,  the  zero  designa¬ 
tion  simply  means  that  the  t*  involving  these  surface  reflection  points  is  too  low  to  be  meas¬ 
ured  by  this  method. 

In  Table  9  and  in  Figure  23a,  the  tp  values  have  been  adjusted  to  reflect  the  t*  values 
expected  for  a  tectonic-tectonic  path  instead  of  the  shield-tectonic  padi  obtained  when  the 
source  event  is  at  East  Kazakh.  A  tp  of  0.15  sec  was  assumed  for  a  shield-shield  path  (Bache 
et  al. ,  1985),  so  half  of  this  was  subtracted  from  die  tp  estimates  for  the  shield-tectonic  padis, 
and  the  remainder  was  doubled,  giving  an  estimate  for  a  tectonic-tectonic  path  where  bodi  tec- 
.  tonic  regions  are  similar  to  the  station  site.  Looking  at  the  tectonic-tectonic  padi,  values  for  tp 

allows  more  direct  comparison  with  the  estimates  of  tpp,  and  the  values  of  tp  and  tpP  are  rela¬ 
tively  close,  though  diey  do  reflect  a  good  deal  of  variability  within  the  tectonic  regions. 

Ihe  PP/P  ratio  method  is,  of  course,  subject  to  errors  due  to  the  distortion  of  the  spectral 
ratios  by  the  effects  ol  crustal  structure  at  the  reflection  points.  We  also  attempted  to  compute 
spectral  ratios  of  PP  to  P  directly  using  the  approach  followed  by  Shore  (1983).  We  have 
lound  that  for  most  events  the  band-pass  approach  was  preferable  because  of  the  poor  S/N 
ratios  of  PP  (noise  means  scattered  P  coda  in  this  case). 

i 

Rise  l  imes  of  P  Waves  from  Nuclear  Explosions 
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Tabic  9 


t*p  and  t*P|i  Estimates  for 
Eurasian  Stations  for  Tectonic  -  Tectonic  Paths 


station 

P  frequency  range 

t*r 

PP  frequency  range 

t*Pp 

(Hz) 

(sec) 

(Hz) 

(sec) 

ANTO 

0.5-4 

0.57±0.08 

0.5-1. 3 

0.5010.09 

CHTO 

0.5-5 

0.43±0.12 

0.5-1. 3 

0.5210.08 

GRFO 

0.5-4 

0.54±0.12 

KAAO 

0.5-6 

0.13±0.05 

KONO 

0.5-5 

0.3410.10 

MAIO 

0.5-4 

0.4110.08 

MAJO 

0.8-5 

0.2510.10 

0.5-1. 3 

0.6010.11 

SHIO 

0.5-5 

0.2910.08 

0.5-1.3 

0.3210.13 

TATO 

1.0-3 

0.1910.09 

52 


Figure  23a.  Map  of  t*  values  for  direct  paths  to  SRO  stations  and  PP  arrivals  (plotted  at 
the  bounce  point).  The  values  of  0  near  two  bounce  points  indicate  low  attenuation  under  the 
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A  time  domain  approach  to  the  estimation  of  Q  for  body  waves  is  the  matching  of  rise 
times  (Stewart,  1984).  This  approach  is  based  on  the  idea  that  t*  is  the  controlling  factor  that 
determines  the  gross  overall  shape  of  the  initial  swing  of  the  short  period  P  waveform  as 
characterized  by  a  "rise  time".  Other  factors  that  may  influence  the  rise  time  are  the  shape  of 
the  initial  source  pulse,  the  time  lag  between  P  and  pP,  and  "stochastic  dispersion"  (Richards 
and  Menke,  1983;  McLaughlin  et  at. ,  1986).  Thus,  by  assuming  that  the  rise  time  is  due 
solely  to  intrinsic  attenuation,  rise  time  estimates  give  upper  bounds  on  t*.  For  nuclear  explo¬ 
sions  we  have  a  reasonably  good  idea  about  the  minimum  rise  times  of  the  source  pulses  from 
near  field  measurements;  these  are  the  shortest  for  shots  in  hard  rock.  The  minimum  P-pP 
time  lag  may  be  obtained  by  assuming  that  the  explosion  was  buried  at  a  relati/ely  shallow 
scaled  depth  (  h=0.07Y1/3,  where  the  depth  h  is  given  in  km  and  the  explosion  yield  Y  in  kilo- 
tons  )  and  assuming  that  the  near  surface  velocities  are  high.  Deeper  burial  depths  and  lower 
uphole  velocities  would  give  overly  high,  conservative  upper  limits  for  tp. 

Rise  time  is  usually  defined  as  the  time  between  the  first  maximum  and  the  intersection 
of  the  tangent  to  the  rising  portion  of  the  P  waveform  with  the  maximum  slope  with  the  zero 
amplitude  level.  Since  the  maximum  slope  is  hard  to  measure  we  have  used  the  time  between 
the  first  break  and  the  first  maximum  as  rise  time.  We  have  used  the  same  definition  in  our 
theoretical  simulations  and,  therefore,  this  modification  of  the  procedure  has  no  effect  on  the 
conclusions.  Theoretical  rise  times  were  computed  using  the  cube-root  scaled  von  Seggern  and 
Blandford  (1972)  granite  source  model,  the  minimum  scale  depth,  a  range  of  apparent  tp  from 
0  to  1  sec,  a  surface  reflection  coefficient  of  negative  unity  and  an  uphole  velocity  of  5.5 
km/sec. 

Rise  times  provide  powerful  arguments  (Stewart,  1984)  against  high  values  of  t  com¬ 
monly  accepted  just  a  few  years  ago  as  valid  for  the  short  period  band.  These  high  values  are 
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not  believed  by  most  seismologists  today  and  the  rise  times  can  only  be  used  to  put  some  loose 
upper  bounds  on  t  .  These  upper  bounds  are  not  very  useful  in  determining  what  the  most 
appropriate  values  of  t  are,  although  they  convincingly  rule  out  the  previously  accepted  high 
values.  We  therefore  relegated  our  rise  time  analyses  to  Appendix  C  and  did  not  repeat  the 
same  kinds  of  analyses  for  the  other  regions  studied  in  this  report.  Nevertheless  the  results  in 
Appendix  C  should  convince  the  reader  that  higher  t  values  are  appropriate  for  "tectonic"  than 
for  "shield"  Eurasia. 

In  Appendix  C,  the  observed  rise  times  are  plotted  against  the  apparent  t*  values  16  for 
some  of  the  SRO  stations  we  used  in  this  study.  We  omitted  KAAO  and  MAIO  because  these 
are  at  regional  distances  from  the  Kazakh  test  site.  The  four  curves  in  each  figure  correspond 
to  the  theoretical  values  for  four  yields  derived  by  measuring  synthetic  seismograms  in  a  simi¬ 
lar  manner.  The  yields  used  to  determine  the  position  of  our  ooservational  data  in  these  graphs 
were  obtained  by  taking  the  mb  values  from  the  seismic  bulletin  published  by  the  International 
Seismic  Centre  and  estimating  the  yields  based  on  the  yield-magnitude  curve  obtained  by  plot¬ 
ting  these  magnitudes  against  the  yields  estimated  by  Dahlman  and  Israelson  (1977)  for 
Kazakh  explosions.  These  yield  estimates,  although  crude,  are  sufficient  for  predicting  the 
pulse  rise  times  since  they  are  very  insensitive  to  yield. 

Each  observed  rise  time  is  plotted  such  that  the  ordinate  is  equal  to  the  observed  rise 
time,  and  the  point  is  located  between  the  curves  to  correspond  to  its  mb-estimated  yield.  The 
projection  of  each  point  on  the  horizontal  axis  gives  an  estimate  of  tp  which,  in  view  of  the 
conservative  assumptions  we  have  used,  must  be  an  upper  limit.  The  results  in  Appendix  C 
indicates  that  the  upper  limits  of  t:!  are  quite  variable  from  station  to  station.  The  lowest  limits 
are  given  by  SHIO  and  ANTO,  intermediate  limits  characterize  KONO  and  GRFO,  while  high 
t  is  associated  with  CHTO.  The  SHIO  results  agree  with  those  of  Stewart  (1984)  who 
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obtained  a  similar  mean  value  for  USSR  explosions  at  YKA,  another  station  situated  on  a 
shield  region. 


The  estimate  tp  from  rise  times  is  significantly  greater  than  the  estimate  tp  generally 
obtained  from  spectra.  Theoretically,  for  a  minimum  phase,  causal  waveform  propagating  in 
an  intrinsically  anelastic  medium,  higher  frequencies  will  travel  faster  than  lower  frequencies. 
This  will  make  the  rise  time  less  than  one  fourth  the  full  cycle  of  the  first  arrival.  Observed 
rise  times,  however,  are  often  equal  to  or  greater  than  one  fourth  the  period  of  the  first  full 
cycle  of  the  P  arrival,  and  it  has  been  observed  that  often  the  higher  frequencies  in  the  P  wave 
arrive  somewhat  later  titan  the  lower  frequencies  (McLaughlin  and  Anderson,  1986;  Richards 
and  Menke,  1983).  This  is  due  to  scattering  which  can  cause  reverse  apparent  dispersion,  with 
the  higher  frequencies  delayed  relative  to  the  lower  frequencies.  Thus,  the  longer  rise  times 
overestimate  t  and  serve  as  upper  bounds  on  the  attenuation.  Recent  theoretical  work  on  sto¬ 


chastic  dispersion  also  indicates  that  such  upper  limits  for  tp  from  rise  times  are  probably 
overestimated  (Richards  and  Menke,  1983;  Frankel  and  Clayton,  1984;  McLaughlin  and  Ander¬ 
son,  1986;  McLaughlin  et  al. ,  1986),  and  therefore  impose  even  more  severe  constraints  on 


the  admissible  maximum  values  of  t  than  previously  thought,  since  the  group  delay  due  to  the 
randomness  of  media  increases  with  increasing  frequency,  thus  increasing  the  rise  times 
independently  of  Q  effects. 

Analyses  of  SS/S  Amplitude  Ratios 

A  set  of  S  and  SS  waveforms  were  collected  at  various  distances.  The  event  parameters 
are  listed  in  Table  10  and  their  locations  and  the  associated  SS  surface  reflection  bounces  are 
shown  in  Figure  23b.  These  have  been  added  to  our  data  for  the  Eurasian  shield  using  the 
approach  in  our  previous  studies  (Der  et  al. ,  1985a).  The  plot  of  our  tectonic  data  points 
(stars)  together  with  our  previous  data  for  the  Eurasian  shield  (filled  circles)  is  shown  in  Figure 
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Table  10 


Parameters  of  earthquakes  used  in  SS/S  study  of  tectonic  Eurasia 


Event  # 
1 
2 

3 

4 


date 

latitude 

longitude 

depth 

mb 

Feb  1,  1984 

49.06 

+  146.59 

573.0 

5.9 

Apr  20,  1984 

50.12 

+  148.75 

582.0 

6.0 

Apr  23,  1984 

47.45 

+  146.69 

414.0 

6.0 

Apr  3,  1985 

28.23 

+139.52 

469.0 

5.9 

I 


riJRASIA 


Figure  23b.  Map  of  the  bounce  points  (+’s)  within  the  tectonic  Eurasia  used  for  the  SS/S 
amplitude  ratio  study. 
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24.  Very  preliminary  results  suggest  high  t*  >  6  sec  for  "tectonic"  Eurasia,  much  higher  than 
the  value  obtained  ('3  sec)  previously  for  the  shield  (Der  et  al. ,  1985b). 

Frequency  Dependent  t*  Model  for  Eurasia 

We  have  obtained  values  of  t  for  both  tectonic  and  shield  Eurasia  using  the  various  data 
sets  and  methods  at  various  different  frequency  bands.  The  short  period  P  wave  studies  pro¬ 
vide  information  towards  the  higher  frequency  band  from  0.6  to  10  Hz.  The  low  frequency 
studies  yield  t  information  between  the  frequency  band  of  .04  to  .1  Hz.  These  results  are 
summarized  in  Figure  25. 

In  constructing  t*(f)  models  for  Eurasia,  we  are  using  the  shield  curve  at  the  higher  fre¬ 
quency  band  as  a  basis  of  reference  in  order  to  estimate  the  Q  structure  under  the  tectonic 
regions  of  Eurasia.  In  Figure  25  we  show  our  t+(f)  curve  for  the  Eurasian  shield  with  some  of 
the  measurements  superposed  for  paths  also  involving  the  tectonically  active  parts  of  Eurasia. 
The  new  data  include  some  of  the  direct  t*  estimates  for  the  WWSSN  stations  located  in  "tec¬ 
tonic"  areas.  Since  these  were  derived  from  Kazakh  explosion  data  the  associated  path  is  a 
shield-tectonic  one.  In  order  to  estimate  what  the  appropriate  t+  for  a  tectonic-tectonic  path 
would  be  we  subtracted  the  shield  contribution  (estimated  at  0.07  sec)  from  these  values  and 
multiplied  the  result  by  two.  These  values  are  plotted  in  the  figure.  The  values  from  PP/P 
ratios  are  directly  applicable  without  adjustments. 

The  results  obtained  in  our  previous  studies  of  the  shield  regions  of  Eurasia  are  plotted 
together  with  the  new  data  for  the  "tectonic"  regions.  The  latter  seems  to  indicate  uniformly 
higher  tp  values  than  for  the  shield.  The  models  are  frequency  dependent  and  are  quasi- 
parallel  relative  to  each  other.  The  total  picture  that  emerges  from  all  this  is  that  the  relative 
attenuation  changes  are  similar  over  the  whole  frequency  band  studied  from  .02  to  10  Hz. 
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ANALYSES  OF  DATA  FROM  NORTH  AMERICA 


Review  of  Previous  Data  Analyses  in  the  Short  and  Intermediate  Period  Bands 


The  attenuation -related  magnitude  bias  has  been  well  studied  in  North  America  and  there 
is  thus  hardly  any  need  to  redo  any  of  the  thorough  analyses  already  in  print.  The  topic  has 
been  discussed  in  a  number  of  papers  (Lay  and  Helmberger,  1981;  Der  era/.,  1980,1982a; 
Der  and  Lees,  1985).  The  spectral  contents  and  amplitudes  of  short  period  P  and  S  waves 
exhibit  some  well  studied  regional  patterns  (Booth  ctal.,  1974;  Lay  and  Helmberger,  1981; 
Der  ct  al.,  1982a).  The  P  wave  amplitudes  have  amplitude  variations  of  up  to  factors  of  2-3 
and  short  period  S  waves  by  factors  of  up  to  8  with  corresponding  large  variations  in  the  dom¬ 
inant  periods  of  these  waves  for  common  events.  It  was  pointed  out  by  Der  and  Lees  (1985) 
that  the  observed  patterns  of  amplitude  and  period  variations  across  North  America  cannot  be 
explained  by  the  models  in  which  the  regional  t*(f)  variations  are  solely  due  to  the  "corner  fre¬ 
quency  shift",  that  is,  by  the  shift  of  the  high  frequency  cutoff  frequency  of  a  single  absorption 
band  model  of  the  type  described  by  Minster  (1978).  Instead,  the  data  require  t*(f)  models  that 
are  quasi-parallel.  In  our  subsequent  syntheses  of  acceptable  models  for  North  America  we 
shall  reconform  this  quasi-parallelism  of  t*(f). 


Summary  of  Results  for  North  America 


Given  the  fact  that  the  relative  mantle  attenuation  is  fairly  well  known  between  the  tec¬ 
tonic  and  shield  areas  of  North  America  in  the  short  and  intermediate  period  bands,  there  is  no 
point  in  redoing  much  of  this  work  for  this  region  (Solomon  and  Toksoz,  1970;  Lay  and  Helm¬ 
berger,  1981;  Solomon,  1972;  Lee  and  Solomon,  1975,1979;  Der  et  al. ,  1982a;  Shore,  1983), 
the  most  important  problem  to  solve  in  this  region  is  for  the  value  of  t  at  low  frequencies. 
These  can  be  obtained  from  the  available  ScS  data.  High  frequency  t*  information  will  he 
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acquired  from  the  numerous  published  data.  Attenuation  estimates  for  die  shield  regions  of 
Noah  America  were  obtained  from  PP/P  spectral  ratios  by  Shore  (1983).  By  studying  15 
events  in  Eurasia  and  recorded  at  N.  American  stations,  Shore  (1983)  computed  spectral  ratios 
of  PP/P  to  arrive  at  an  average  tp  of  approximately  0.29  sec  for  shield  N.  America.  In  addi¬ 
tion,  Der  et  al.  (1982a)  have  also  obtained  tp  in  the  range  of  0.1  to  0.2  sec  for  shield  paths 
and  tp  differential  between  shield-to-shield  and  shield-to-tectonic  paths  to  be  of  the  order  of 
0.2  s.  We  shall  adopt  these  results  for  the  shoa  period  band  to  integrate  with  those  obtained 
for  the  low  frequencies  using  ScS  phases  to  construct  frequencies  dependent  t*  curves  for  N. 
America. 

By  combining  t*  values  covering  the  short  and  long  period  bands,  we  construct  t*(f) 
curves  for  tectonic  and  shield  N.  America  and  are  plotted  together  in  Figure  26.  In  this  figure 
we  show  our  best  estimates  of  the  t*(f)  curves  that  fit  our  short  and  long  period  analyses.  In 
compiling  this  figure  we  have  doubled  die  SWUS-shield  t  differentials  to  obtain  the  appropri¬ 
ate  tectonic-tectonic  path  t  estimates.  While  the  data  for  the  southwestern  US  show  a  clear 
frequency  dependence  of  t\  frequency  dependence  cannot  be  shown  for  the  Canadian  shield 
from  these  data.  These  curves  indicate  the  same  quasi-parallel  pattern  as  that  observed  for 
Eurasia,  with  the  tectonic  region  being  higher  than  the  shield  region.  The  value  of  t*  decreases 
as  the  frequencies  increases  in  all  cases.  It  is  obvious  that  the  two  curves  are  well  separated 
and  any  systematic  bias  will  not  affect  their  relative  shapes. 
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ANALYSES  OF  DATA  FROM  EASTERN  PACIFIC 

Review  of  Previous  Data  Analyses  in  tiie  Short  and  Intermediate  Period  Bands 

Large  sets  of  P  wave  data  have  been  analyzed  from  French  nuclear  explosions  from 
Polynesia  by  several  workers  (Der  et  al. ,  1985b;  Bache  et  al. ,  1985).  The  consensus  is  that 
P  waves  observed  at  teleseismic  distances  from  this  region  have  spectral  characteristics  similar 
to  those  from  NFS  thus  indicating  a  low  Q  under  the  Tuamoto  Archipelago.  The  t*  contribu¬ 
tion  to  a  teleseismic  P  wave  by  the  upper  mantle  under  Tuamoto  was  estimated  as  0.35  sec  by 
Der  et  al. ,  (1985b)  and  by  McLaughlin  et  al. ,  (1986)  to  be  around  0.43  sec. 

PP/PKP  Spectral  Ratios 

To  study  the  short  period  attenuation  for  the  S.  Pacific,  we  employ  the  similar  procedure 
as  used  in  an  earlier  section.  Shapes  of  spectral  ratio  of  PP-PKP  pairs  are  computed  to  obtain 
tp  across  the  S.  Pacific  Ocean.  The  PP  phase  would  bounce  off  the  surface  and  sample  the 
upper  mantle  near  the  bounce  point  with  two  passages  whereas  PKP  phase  would  diffract  off 
the  core  and  contains  no  reflection  from  the  surface  of  the  Earth.  The  ratio  between  PP-PKP 
would  provide  an  estimate  of  upper  mantle  attenuation  since  the  Q  in  the  lower  mantle  and 
core  is  high  compared  to  that  of  the  upper  mantle.  These  phases  are  bandpass  filtered  in  nar¬ 
row  frequency  bands  between  0.2  to  6  Hz.  The  bandpass  filters  allow  rapid  identification  of 
the  phases  with  estimates  of  the  noise  level.  Amplitudes  are  measured  in  the  time  domain  for 
each  phase  and  PP/PKP  amplitude  ratios  are  computed  for  each  frequency  band. 

Bi-directional  trans-S.  Pacific  paths  are  collected  for  deep  focus  Fiji-Tonga  earthquakes 
recoided  at  S.  American  stations  and  for  deep  focus  S.  American  earthquakes  recorded  at 
southwestern  Pacific  stations.  We  have  examined  mainly  deep  focus  events  associated  with  the 
subducting  slabs  in  S.  America  and  Tonga-Fiji  regions  recorded  using  digital  DWWSSN  and 
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SRO  stations.  In  most  of  these  events,  the  PP  are  either  highly  attenuated  or  are  contaminated 
by  other  arrivals. 

Short  period  vertical  component  recordings  from  twelve  deep  focus  earthquakes  are 
assembled  to  study  tp  for  the  S.  Pacific  Ocean  using  PP/PKP  amplitude  ratio  method.  The 
locations  of  these  events  and  stations  are  indicated  in  Figure  27.  The  event  parameters  are 
listed  in  Table  11.  The  surface  bounce  points  for  the  PP  phases  are  also  plotted  in  the  same 
figure.  All  of  the  PP  phases  have  surface  reflection  bounce  points  within  die  S.  Pacific  and  are 
well  away  from  the  axis  of  the  E.  Pacific  Rise.  Due  to  the  long  epicentral  distance  across  the 
Pacific,  only  multiple  PP  and  PKP  phases  are  observable  as  all  P  phases  are  diffracted.  Since 
the  earthquakes  used  are  deep,  the  ratio  of  PP  and  PKP  approximately  represents  die  anelasti- 
city  properties  of  the  upper  mantle  in  canceling  out  the  source  effects. 

A  similar  set  of  band-pass  filters  are  applied  to  our  data  set  here  and  the  filtered  plots  are 
shown  in  Appendix  B  (Figures  B9-B20).  The  arrival  times  for  both  phases  are  identified  using 
the  J-B  travel  time  table.  The  frequency  band  of  PKP  ranges  from  0.3  to  3  Hz  while  PP  phase 
generally  emerges  at  a  lower  frequency  ranging  from  0.1  to  2.5  Hz.  This  form  of  display 
allows  one  to  measure  the  amplitudes  within  good  signal  to  noise  level. 

The  amplitude  of  these  phases  is  measured  from  a  frequency  band  of  0.4  to  6  Hz.  Where 
there  is  too  low  a  signal  in  either  phase,  no  measurements  are  attempted  to  constraint  the  curve 
giving  rise  to  bias  in  constructing  the  curve.  Amplitude  ratios  for  each  event  are  plotted  in 
linear-log  scale  and  an  example  is  shown  in  Figure  28.  A  least  square  fit  to  the  data  points 
gives  a  best  fitting  curve  through  the  PP/PKP  amplitude  ratio  data.  By  averaging  the  tpp  for 
all  events,  we  obtained  a  tpp  of  0.84  ±  0.30  sec  for  die  S.  Pacific  in  close  agreement  to  the 
value  obtained  by  Der  ct  al.,  (1985b)  and  McLaughlin  ct  al.,  (1986). 
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Parameters  of  earthquakes  used  in  PP/PKP  study  of  tectonic  Pacific 


Event  # 

date 

1 

Nov  25,  1976 

2 

Sep  6,  1978 

3 

Sep  23,  1978 

4 

Nov  8,  1979 

5 

Jan  4,  1982 

6 

Sep  16,  1983 

7 

Nov  29,  1983 

8 

Apr  22,  1984 

9 

May  30,  1984 

10 

Jun  15,  1984 

11 

Nov  17,  1984 

12 

Oct  6,  1985 

latitude 

longitude 

depth 

mb 

-19.50 

-177.70 

424 

6.0 

-13.35 

+  167.13 

214 

5.9 

-13.94 

+167.18 

216 

6.1 

-32.26 

+179.33 

447 

5.7 

-23.17 

-177.31 

195 

6.0 

-24.03 

-179.80 

510 

6.0 

-19.50 

-177.78 

525 

5.7 

-21.87 

-179.38 

593 

5.7 

-04.84 

+  151.58 

174 

6.2 

-15.82 

-174.83 

247 

6.1 

-18.79 

-178.03 

451 

6.1 

-18.96 

169.43 

273 

5.7 

68 
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SSS/SS  Amplitude  Ratios 

In  order  to  study  the  variations  of  long  period  t  in  Eastern  Pacific,  we  have  employed  an 
amplitude  ratio  method  between  SS  and  SSS  phases.  The  same  approach  has  been  used  in 
other  regions  as  discussed  in  our  previous  studies  (Der  et  al.,  1985a).  The  SS  phase  has  one 
surface  bounce  point  whereas  the  SSS  phase  has  two  surface  bounces  along  the  path.  The 
amplitude  ratio  between  the  two  phases  indicates  the  energy  loss  caused  by  a  double  passage 
of  the  S  wave  through  the  upper  mande  anelasdc  structure. 

We  have  examined  deep  focus  events  in  Fiji-Tonga  and  S.  America  recorded  at  S.  Ameri¬ 
can  and  Western  Pacific  stations  respectively.  Both  the  SS  and  SSS  are  very  difficult  to  iden¬ 
tify  due  to  the  relatively  low  Q  in  the  tectonic  Pacific  and  contaminations  due  to  other  arrivals. 
The  locations  of  the  earthquake  source  and  station  used  in  the  SS-SSS  measurement  are  shown 
in  Figure  29  and  the  parameters  are  listed  in  Table  12.  The  location  of  the  bounce  points 
shown  in  the  figure  indicate  that  they  do  not  sample  the  East  Pacific  Rise  spreading  center  but 
are  scattered  among  various  age  regions  within  the  Pacific.  Since  our  interest  is  to  estimate  t* 
for  the  broad  E.  Pacific,  at  this  point  we  do  not  attempt  to  resolve  the  variation  of  Q  within  the 
different  age  regions  due  to  the  scatter  of  our  data. 

The  amplitude  ratios  of  the  SSS/SS  are  plotted  in  Figure  30  fitted  with  various  theoretical 
ts  curves.  These  theoretical  curves  are  computed  assuming  an  oceanic  model  and  a  tg  of  6  sec 
for  the  SS  path.  The  5t  between  SS  and  SSS  varies  from  1  to  7  sec.  From  a  rms  fit  to  the 
data,  we  arrive  at  an  estimate  of  6  sec  for  of  shear  waves  across  the  E.  Pacific.  The  values 
obtained  here  provide  a  constraint  to  the  long  period  data  in  constructing  a  t*  curve  for  the 
Pacific. 
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Parameters  of  earthquakes  used  in  SSS/SS  study  of  tectonic  Pacific 


Event  # 

date 

latitude 

longitude 

depth 

mb 

1 

Oct  22,  1977 

-27.96 

-62.82 

592.3 

6.0 

2 

May  13,  1978 

-14.60 

+  167.31 

181.7 

6.5 

3 

Jul  17,  1978 

-14.95 

-175.81 

305.4 

5.9 

4 

Jul  17,  1978 

-14.95 

-175.81 

305.4 

5.9 

5 

Sep  6,  1978 

-13.35 

+  167.13 

214.8 

5.9 

6 

Sep  23,  1978 

-13.94 

+  167.18 

216.7 

6.1 

7 

Apr  24,  1979 

-20.82 

-178.67 

596.1 

5.9 

8 

Aug  5,  1979 

-22.86 

-177.47 

203.7 

6.1 

9 

Apr  13,  1980 

-23.53 

-177.27 

102.1 

6.6 

10 

Jul  20,  1980 

-17.88 

-178.61 

587.8 

6.0 

11 

Nov  30,  1980 

-19.43 

-175.85 

202.0 

6.0 

12 

Nov  30,  1980 

-19.43 

-175.85 

202.0 

6.0 

13 

Sep  28,  1981 

-29.33 

-179.04 

323.0 

6.0 

14 

Nov  25,  1981 

-20.94 

-178.86 

614.0 

5.9 

15 

Sep  16,  1983 

-24.03 

-179.80 

510.0 

6.0 

'6 

Dec  21,  1983 

-28.19 

-63.17 

602.0 

6.2 

17 

May  1,  1985 

-09.20 

-71.23 

600.0 

6.0 
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Final  Model  for  E.  Pacific 

In  general,  the  broad  Pacific  Ocean  is  characterized  as  tectonic,  therefore  only  one  curve 

* 

will  be  constructed  for  the  region,  t  values  for  the  high  frequencies  are  obtained  from  PP/PKP 
spectral  ratios  whereas  the  low  frequency  information  is  obtained  from  multiply  reflected  ScS 
phases  and  SS/SSS  amplitude  ratios.  Combining  the  three  sets  of  data,  a  broad-band  t*(f) 
model  is  constructed  for  the  E.  Pacific  as  shown  in  Figure  31.  This  model  is  frequency  depen¬ 
dent  with  higher  t+  at  the  lower  frequencies  of  0.02  Hz  and  decreases  towards  the  higher  fre¬ 
quencies  of  up  to  2  Hz,  similar  in  form  to  the  ones  obtained  earlier  for  the  other  regions.  The 
long  period  estimates  front  ScS  and  SS/SSS  are  not  in  exact  agreement  to  each  other.  This  is 
most  likely  due  to  the  scatter  of  bounce  points  within  the  E.  Pacific  and  the  different  frequency 
contents  of  these  phases.  In  constructing  the  curves,  more  emphasis  is  being  put  on  using  the 
ScSN  data  since  they  are  generally  of  higher  quality  than  the  SS  and  SSS  data. 
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DISCUSSION  OF  THE  BROAD  BAND  RESULTS 

Reviewing  the  ScS  results  together  with  die  data  presented  for  the  individual  regions  of 
study  earlier  in  this  report,  the  pattern  is  evident:  wherever  the  anelastic  attenuation  is  shown  to 
be  high  in  the  short  period  and  intermediate  period  bands,  it  is  also  high  in  the  long  period 
band  spanned  by  multiple  S  and  ScS  data.  In  areas  where  the  tp  is  low,  the  ScS  phases  have 
also  a  lower  attenuation,  the  PP/P  have  a  higher  percentage  of  high  frequencies  and  the  rise 
times  of  P  and  S  are  also  shorter.  In  areas  of  high  attenuation  just  the  reverse  is  true.  The 
variations  in  wave  characteristics  in  our  "shield"  and  "tectonic"  regions  are  similar  in  nature, 
but  not  as  severe,  to  those  observed  in  areas  of  anomalously  high  and  low  attenuation,  such  as 
spreading  back-arc  basins  and  spreading  ridge  axis  respectively.  This  way  our  results  fit  into 
the  framework  of  the  picture  provided  by  the  "new  global  tectonics".  In  a  way  this  is  to  be 
expected,  since  regional  variations  of  Q  are  thought  to  be  caused  in  a  large  part  by  lateral  vari¬ 
ations  of  temperature  thus  causing  variations  of  Q  which  are  of  similar  sign  over  wide  fre¬ 
quency  bands. 

We  summarize  in  Table  13  the  t*  or  t*  estimates,  whiche  er  one  is  the  most  applicable, 
for  the  various  types  of  measurements  and  frequency  bands  in  all  the  regions  studied.  These 
results  give  the  raw  material  for  inversions  for  Q(f,r)  and  the  definition  of  the  "best"  t  (0 
models. 
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Table  13 

Comparison  of  long  period  ScS  and  short  period  P 
attenuation  results 


Region 

ts*cS  (.02-.  1  Hz) 

tp(  .5-10  Hz) 

Tectonic  Eurasia 

3.0-3. 9 

0.5 

Shield  Eurasia 

2.5 

0.15-0.2 

Tectonic  N.  America 

4.0-4.5 

0.4 

Shield  N.  America 

1. 8-2.5 

0.15 

SE.  tran-Pacific 

5.5-7.0 

0.8 

W.  Pacific 

4.0-5.0 

Shield  S.  America 


2. 0-3.0 
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POSSIBLE  t*(0  and  Q(f,r)  MODELS  AND  YIELD  ESTIMATION 

We  are  now  ready  to  outline  the  t*(f)  curves  most  appropriate  to  each  of  the  regions  stu¬ 
died.  In  Figure  32  we  show  these  curves  that  have  been  obtained  based  on  the  data  in  Table 
13.  These  curves  apply  to  teleseismic  paths  confined  entirely  to  the  respective  structures.  In 
arriving  at  a  gross  estimates  of  t*(f),  we  do  not  attempt  to  resolve  the  lateral  variation  of  Q 
within  each  region,  for  example,  the  different  age  regions  within  the  E.  Pacific.  Curves  for 
mixed  paths  can  be  constructed  by  linearly  combining  pairs  of  curves.  The  results  do  not 
favor  the  types  of  convergent  or  divergent  models  of  t*(f)  that  could  confound  the  spectral 
method  in  estimating  relative  attenuation  corrections  derived  from  spectra  (Hadley  and  Mell- 
man,  1983;  Lay  and  Helmberger,  1981),  rather  they  favor  quasi-parallelism  (Der  and  Lees, 
1985).  It  is  quite  possible  that  more  subtle,  second  order  errors  in  attenuation  may  occur  by 
weak  convergence  or  divergence  of  such  frequency  dependent  t*  curves.  We  also  neglect  here 
any  distance  dependence  of  t\  yet  another  second  order  effect. 

Our  data  do  not  allow  us  to  make  a  strict  inversion  for  Q(f,r),  since  the  depth  dependence 
is  not  well  constrained.  Nevertheless,  assuming  that  the  lateral  Q  variations  in  the  Earth  is  the 
most  pronounced  within  a  200-300  km  thick  layer  in  the  upper  mantle,  we  can  confidently  state 
that  Q  may  vary  by  about  a  factor  of  5-10  laterally  over  wide  frequency  range  of  .01  to  .1  Hz. 
Such  a  strong  lateral  variation  in  Q  indicates  that  it  is  highly  unlikely  that  t  (0  models  con¬ 
verge  at  low  frequencies. 

With  regards  to  yield  estimation,  the  postulated  types  of  t’(f)  allow  one  to  predict  relative 
attenuation  corrections  for  various  source  and  receiver  regions  derived  from  shapes  of  P  and  S 
wave  spectra  as  done  before  (Der  ct  al. ,  1985b).  Mixed  path  attenuation  may  also  be  corrected 
by  using  a  linear  combination  of  the  t*(f)  values  for  different  regions.  On  the  other  hand, 
since  our  data  require  that  the  Q  be  frequency  dependent  in  all  the  regions  studied  the  absolute 
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signal  (P  and  S  wave)  amplitudes  will  be  altered  by  the  frequency-dependent  Q  models  relative 
to  those  predicted  by  constant  Q  derived  from  short  period  data.  This  removes  at  least  some 
of  the  apparent  (but  not  real)  discrepancies  between  absolute  amplitudes  and  the  so  called 
"spectral  Q"  claimed  (Burdick  ct  al. ,  1984). 
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CONCLUSIONS 

Broad-band  models  of  t*(f)  are  presented  for  different  tectonic  and  shield  regions.  These 
regions  include  tectonic  N.  America,  Eurasia,  E.  Pacific,  W.  Pacific  and  S.  America  and  shield 
N.  America,  Eurasia,  and  S.  America.  These  models  are  discussed  in  conjunction  with  each 
other  in  terms  of  the  lateral  variation  and  frequency  dependence  of  Q  in  various  parts  of  the 
world.  The  main  question,  from  the  point  of  view  of  practical  consequences  of  this  study  in 
estimating  attenuation  is  the  approximate  form  of  t  (f)  for  teleseismic  paths  crossing  the  mantle 
under  various  regions. 

In  Eurasia  we  are  using  the  shield  curve  obtained  previously  as  a  basis  of  reference  in 
order  to  estimate  the  Q  structure  under  the  tectonic  regions  of  Eurasia.  We  have  also  presented 
some  direct  t  estimates  for  the  WWSSN  stations  located  in  "tectonic"  areas.  All  of  the  obser¬ 
vations  for  the  tectonic  areas  are  higher  than  the  corresponding  measurements  under  shields. 
The  total  picture  that  emerges  from  all  this  is  that  the  relative  attenuation  changes  among  the 
two  regions  are  similar  in  sign  over  the  whole  frequency  band  studied  from  .02  to  10  Hz. 

T  here  is  no  question  that  in  the  short  period  band  at  frequencies  above  2  Hz,  Q  is  quite 
high  under  the  shield  regions  of  Eurasia  (Der  and  McElfrcsh,  1977;  Der  et  al. ,  1982b)  with 
values  not  less  than  500  for  Qp.  The  argument  of  Der  and  Lees  (1985)  uses  the  assumption 
of  =4  tp.  Their  conclusion  that  the  two  t  (f)  curves  do  not  converge  at  low  frequencies 
would  even  be  strengthened  if  we  assume  a  lower  ratio  between  the  t*  appropriate  to  S  waves 
and  P  waves  (Taylor  et  al.,  1985).  We  do  not  believe,  however,  that  the  ratio  is  considerably 
less  than  4. 

A  considerable  amount  of  work  has  already  been  done  on  studying  attenuation  under  con¬ 
tinental  North  America.  In  this  study  we  compile  the  results  of  past  work  and  supplement 
tilt,  m  with  additional  long  period  data  and  address  the  issue  in  a  broad-band  context.  The  most 
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controversial  question  concerns  the  changes  in  Q  applicable  to  long  period  band  data.  Solo¬ 
mon  and  Toksoz  (1970)  and  Lee  and  Solomon  (1975,1979)  found  large  regional  variations  in 
Q  in  the  long  period  band,  while  Lay  and  Heimberger  (1981)  expressed  doubts  about  any  such 
variations.  In  the  short  and  intermediate  period  bands  Taylor  et  al.  (1985)  and  Der  and  Lees 
(1985)  have  studied  P  and  S  wave  attenuation  and  found  that  significant  Q  differentials  existed 
in  the  upper  mantle  among  various  regions  of  the  United  States  that  were  detectabh  throughout 
the  total  range  of  the  frequency  bands  utilized  in  these  studies. 

Der  and  Lees  (1985)  investigated  some  t*(f)  curves  that  were  proposed  for  purely  shield 
and  shield-to-tectonic  paths  involving  the  mantle  of  North  America.  Using  spectral  and 
waveform  data,  their  study  showed  that  only  pairs  of  curves  that  are  quasi-parallel  (QP)  can 
explain  the  results  of  both  kinds  of  analyses.  Several  models  that  converged  towards  low  fre¬ 
quencies  can  be  rejected  based  on  this  kind  of  data.  It  was  also  concluded  that  using  the  short 
period  data  alone,  it  is  not  possible  to  specify  the  exact  forms  of  frequency  dependence  since 
most  pairs  of  quasi-parallel  t+(f)  curves  can  explain  the  kind  of  short  period  P  and  S  data 
presented  by  Der  and  Lees  (1985). 

The  E.  Pacific  attenuation  model  has  the  highest  t*  among  all  regions.  The  curve  also 
assumes  the  same  frequency  dependent  pattern  as  the  other  regioas  being  quasi-parallel  to 
them.  The  low  frequency  t  for  the  E.  Pacific  is  up  to  a  factor  of  three  larger  than  that  for  the 
shield  regions,  strongly  indicating  that  t  does  not  converge  towards  the  low  trequencies.  The 
t*(f)  curve  being  higher  than  those  for  tectonic  N.  America  and  Eurasia  suggests  that  regional 
variation  in  attenuation  among  tectonic  regions  are  also  inevitable.  We  have  imposed  con¬ 
straints  on  our  work  by  studying  other  regions  and  wherever  possible  reprocess  data  that  were 
used  in  other  published  studies.  These  regions  include  tectonic  and  shield  S.  America,  W.  and 
NW.  Pacific.  Results  from  our  study  are  similar  to  others,  and  being  higher  in  some  cases. 
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Lateral  variations  in  low  frequency  tscS  are  observed  in  these  additional  regions  in  support  of 
our  earlier  conclusion. 

Since  our  data  do  not  have  the  resolution  to  warrant  a  full  inversion  on  the  t*(f)  models, 
to  obtain  Q(r,f)  we  can  only  say  that  the  finding  of  Solomon  and  Toksoz  (1970)  that  large  Q 
differentials  exist  in  the  upper  mantle  at  low  frequencies  across  North  America  seems  to  be 
supported  by  our  data  thus  far.  The  same  is  true  for  all  other  regions  we  have  studied.  Since 
the  ratios  of  t  are  less  than  the  ratios  in  Q  in  the  upper  mantle  in  order  to  cause  the  observed 
amplitude  decay  and  spectral  differences  (in  proportion  to  the  travel  times  through  the  total 
mantle  and  the  upper  mantle). 

In  addition,  in  all  the  regions  we  have  studied,  the  t*  estimates  tend  to  rise  towards  lower 
frequencies  supporting  the  idea  of  frequency  dependence  of  Q  (Lundquist  and  Cormier,  1978; 
Sipkin  and  Jordan,  1979;  Anderson  and  Given,  1982).  Invariably  in  the  short  period  band, 
more  high  frequency  energy  is  seen  than  that  compatible  with  the  Q  estimates  from  long  period 
data.  With  regards  to  yield  estimation  these  results  are  important  since  they  diminish  the  likeli¬ 
hood  of  the  kinds  of  non-parallel  t  (f)  models  that  would  make  impossible  estimation  of  mag¬ 
nitude  bias  from  spectral  data. 
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APPENDIX  A 


Time  Domain  Analyses  of  Multiple  ScS  Phases. 
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Analyses  of  PP/P  and  PP/PKP  Spectral  Ratios. 
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Figure  Bl.  Band-pass  filtered  P  and  PP  phases  on  vertical  component  for  a  Kazakh  nuclear 
explosion  at  KONO.  The  PP  phase  is  observable  up  to  6  Hz  indicating  high  Q  surface 
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Figure  B3.  Band-pass  filtered  P  and  PP  phases  on  vertical  component  for  a  Kazakh  nuclear 
explosion  at  ANTO.  The  PP  phase  is  observable  up  to  1.8  Hz. 
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Figure  B4.  Band-pass  filtered  P  and  PP  phases  on  vertical  component  for  a  Kazakh  nuclear 
explosion  at  SHIO.  The  PP  phase  is  questionable. 


u~n  vk  v>  v*  vx  w  wx  u*  .  j*  jvw  a  wiv  uva  n"  mw  s"  ^  sn  v  o  v  v\  v  c  k/\  k.^  vi  vi .  v  v  > .  xwv  a  ^  r  n  > ..  ■  o  r« 


Figure  B7.  Band-pass  filtered  P  and  PP  phases  on  vertical  component  for  a  Kazakh  nuclear 
explosion  at  CHTO.  The  PP  phase  may  not  be  seen  above  0.8  Hz. 
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Figure  Bll.  Band-pass  filtered  vertical  component  seismograms  for  path  across  E.  Pacific 
with  arrival  times  of  PKP  and  PP  identified.  Absence  of  high  frequency  PP  phase  indicates 
a  low  Q  bounce  region. 
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Figure  B17.  Band-pass  filtered  vertical  component  seismograms  for  path  across  E.  Pacific 

with  arrival  times  of  PKP  and  PP  identified.  Absence  of  high  frequency  PP  phase  indicates 

a  low  Q  bounce  region. 
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Figure  B18.  Band-pass  filtered  vertical  component  seismograms  for  path  across  E.  Pacific 
with  arrival  times  of  PKP  and  PP  identified.  Absence  of  high  frequency  PP  phase  indicates 
!  a  low  Q  bounce  region. 
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